An immune cell atlas reveals the dynamics of human macrophage specification during prenatal development.

Macrophages are heterogeneous and play critical roles in development and disease, but their diversity, function, and specification remain inadequately understood during human development. We generated a single-cell RNA sequencing map of the dynamics of human macrophage specification from PCW 4-26 across 19 tissues. We identified a microglia-like population and a proangiogenic population in 15 macrophage subtypes. Microglia-like cells, molecularly and morphologically similar to microglia in the CNS, are present in the fetal epidermis, testicle, and heart. They are the major immune population in the early epidermis, exhibit a polarized distribution along the dorsal-lateral-ventral axis, and interact with neural crest cells, modulating their differentiation along the melanocyte lineage. Through spatial and differentiation trajectory analysis, we also showed that proangiogenic macrophages are perivascular across fetal organs and likely yolk-sac-derived as microglia. Our study provides a comprehensive map of the heterogeneity and developmental dynamics of human macrophages and unravels their diverse functions during development.

Paradoxical Mitophagy Regulation by PINK1 and TUFm.

Aberrant mitophagy has been implicated in a broad spectrum of disorders. PINK1, Parkin, and ubiquitin have pivotal roles in priming mitophagy. However, the entire regulatory landscape and the precise control mechanisms of mitophagy remain to be elucidated. Here, we uncover fundamental mitophagy regulation involving PINK1 and a non-canonical role of the mitochondrial Tu translation elongation factor (TUFm). The mitochondrion-cytosol dual-localized TUFm interacts with PINK1 biochemically and genetically, which is an evolutionarily conserved Parkin-independent route toward mitophagy. A PINK1-dependent TUFm phosphoswitch at Ser222 determines conversion from activating to suppressing mitophagy. PINK1 modulates differential translocation of TUFm because p-S222-TUFm is restricted predominantly to the cytosol, where it inhibits mitophagy by impeding Atg5-Atg12 formation. The self-antagonizing feature of PINK1/TUFm is critical for the robustness of mitophagy regulation, achieved by the unique kinetic parameters of p-S222-TUFm, p-S65-ubiquitin, and their common kinase PINK1. Our findings provide new mechanistic insights into mitophagy and mitophagy-associated disorders.

TransIntegrator: capture nearly full protein-coding transcript variants via integrating Illumina and PacBio transcriptomes.

Genes have the ability to produce transcript variants that perform specific cellular functions. However, accurately detecting all transcript variants remains a long-standing challenge, especially when working with poorly annotated genomes or without a known genome. To address this issue, we have developed a new computational method, TransIntegrator, which enables transcriptome-wide detection of novel transcript variants. For this, we determined 10 Illumina sequencing transcriptomes and a PacBio full-length transcriptome for consecutive embryo development stages of amphioxus, a species of great evolutionary importance. Based on the transcriptomes, we employed TransIntegrator to create a comprehensive transcript variant library, namely iTranscriptome. The resulting iTrancriptome contained 91 915 distinct transcript variants, with an average of 2.4 variants per gene. This substantially improved current amphioxus genome annotation by expanding the number of genes from 21 954 to 38 777. Further analysis manifested that the gene expansion was largely ascribed to integration of multiple Illumina datasets instead of involving the PacBio data. Moreover, we demonstrated an example application of TransIntegrator, via generating iTrancriptome, in aiding accurate transcriptome assembly, which significantly outperformed other hybrid methods such as IDP-denovo and Trinity. For user convenience, we have deposited the source codes of TransIntegrator on GitHub as well as a conda package in Anaconda. In summary, this study proposes an affordable but efficient method for reliable transcriptomic research in most species.

SperMD: the expression atlas of sperm maturation.

The impairment of sperm maturation is one of the major pathogenic factors in male subfertility, a serious medical and social problem affecting millions of global couples. Regrettably, the existing research on sperm maturation is slow, limited, and fragmented, largely attributable to the lack of a global molecular view. To fill the data gap, we newly established a database, namely the Sperm Maturation Database (SperMD, http://bio-add.org/SperMD ). SperMD integrates heterogeneous multi-omics data (170 transcriptomes, 91 proteomes, and five human metabolomes) to illustrate the transcriptional, translational, and metabolic manifestations during the entire lifespan of sperm maturation. These data involve almost all crucial scenarios related to sperm maturation, including the tissue components of the epididymal microenvironment, cell constituents of tissues, different pathological states, and so on. To the best of our knowledge, SperMD could be one of the limited repositories that provide focused and comprehensive information on sperm maturation. Easy-to-use web services are also implemented to enhance the experience of data retrieval and molecular comparison between humans and mice. Furthermore, the manuscript illustrates an example application demonstrated to systematically characterize novel gene functions in sperm maturation. Nevertheless, SperMD undertakes the endeavor to integrate the islanding omics data, offering a panoramic molecular view of how the spermatozoa gain full reproductive abilities. It will serve as a valuable resource for the systematic exploration of sperm maturation and for prioritizing the biomarkers and targets for precise diagnosis and therapy of male subfertility.

ROS-induced oxidative stress is a major contributor to sperm cryoinjury.

STUDY QUESTION: What is the mechanism behind cryoinjury in human sperm, particularly concerning the interplay between reactive oxygen species (ROS) and autophagy, and how does it subsequently affect sperm fate? SUMMARY ANSWER: The freeze-thaw operation induces oxidative stress by generating abundant ROS, which impairs sperm motility and activates autophagy, ultimately guiding the sperm toward programmed cell death such as apoptosis and necrosis, as well as triggering premature capacitation. WHAT IS KNOWN ALREADY: Both ROS-induced oxidative stress and autophagy are thought to exert an influence on the quality of frozen-thawed sperm. STUDY DESIGN, SIZE, DURATION: Overall, 84 semen specimens were collected from young healthy fertile males, with careful quality evaluation. The specimens were split into three groups to investigate the ROS-induced cryoinjury: normal control without any treatment, sperm treated with 0.5 mM hydrogen peroxide (H2O2) for 1 h, and sperm thawed following cryopreservation. Samples from 48 individuals underwent computer-assisted human sperm analysis (CASA) to evaluate sperm quality in response to the treatments. Semen samples from three donors were analyzed for changes in the sperm proteome after H2O2 treatment, and another set of samples from three donors were analyzed for changes following the freeze-thaw process. The other 30 samples were used for fluorescence-staining and western blotting. PARTICIPANTS/MATERIALS, SETTING, METHODS: Sperm motility parameters, including progressive motility (PR %) and total motility (PR + NP %), were evaluated using the CASA system on a minimum of 200 spermatozoa. The proteomic profiles were determined with label-free mass spectrometry (MS/MS) and protein identification was performed via ion search against the NCBI human database. Subsequently, comprehensive bioinformatics was applied to detect significant proteomic changes and functional enrichment. Fluorescence-staining and western blot analyses were also conducted to confirm the proteomic changes on selected key proteins. The ROS level was measured using 2',7'-dichlorodihydrofluorescein diacetate labeling and the abundance of bioactive mitochondria was determined by evaluating the inner mitochondrial membrane potential (MMP) level. Molecular behaviors of sequestosome-1 (p62 or SQSTM1) and microtubule-associated proteins 1A/1B light chain 3 (LC3) were monitored to evaluate the state of apoptosis in human sperm. Fluorescent probes oxazole yellow (YO-PRO-1) and propidium iodide (PI) were utilized to monitor programmed cell death, namely apoptosis and necrosis. Additionally, gradient concentrations of antioxidant coenzyme Q10 (CoQ10) were introduced to suppress ROS impacts on sperm. MAIN RESULTS AND THE ROLE OF CHANCE: The CASA analysis revealed a significant decrease in sperm motility for both the H2O2-treatment and freeze-thaw groups. Fluorescence staining showed that high ROS levels were produced in the treated sperm and the MMPs were largely reduced. The introduction of CoQ10 at concentrations of 20 and 30 µM resulted in a significant rescue of progressive motility (P < 0.05). The result suggested that excessive ROS could be the major cause of sperm motility impairment, likely by damaging mitochondrial energy generation. Autophagy was significantly activated in sperm when they were under oxidative stress, as evidenced by the upregulation of p62 and the increased conversion of LC3 as well as the upregulation of several autophagy-related proteins, such as charged multivesicular body protein 2a, mitochondrial import receptor subunit TOM22 homolog, and WD repeat domain phosphoinositide-interacting protein 2. Additionally, fluorescent staining indicated the occurrence of apoptosis and necrosis in both H2O2-treated sperm and post-thaw sperm. The cell death process can be suppressed when CoQ10 is introduced, which consolidates the view that ROS could be the major contributor to sperm cryoinjury. The freeze-thaw process could also initiate sperm premature capacitation, demonstrated by the prominent increase in tyrosine phosphorylated proteins, verified with anti-phosphotyrosine antibody and immunofluorescence assays. The upregulation of capacitation-related proteins, such as hyaluronidase 3 and Folate receptor alpha, supported this finding. LARGE SCALE DATA: The data underlying this article are available in the article and its online supplementary material. LIMITATIONS, REASONS FOR CAUTION: The semen samples were obtained exclusively from young, healthy, and fertile males with progressive motility exceeding 60%, which might overemphasize the positive effects while possibly neglecting the negative impacts of cryoinjury. Additionally, the H2O2 treatment conditions in this study may not precisely mimic the oxidative stress experienced by sperm after thawing from cryopreservation, potentially resulting in the omission of certain molecular alterations. WIDER IMPLICATIONS OF THE FINDINGS: This study provides substantial proteomic data for a comprehensive and deeper understanding of the impact of cryopreservation on sperm quality. It will facilitate the design of optimal protocols for utilizing cryopreserved sperm to improve applications, such as ART, and help resolve various adverse situations caused by chemotherapy, radiotherapy, and surgery. STUDY FUNDING/COMPETING INTEREST(S): This work was supported by grants from the Major Innovation Project of Research Institute of National Health Commission (#2022GJZD01-3) and the National Key R&D Program of China (#2018YFC1003600). All authors declare no competing interests. TRIAL REGISTRATION NUMBER: N/A.

Mining Real-World Big Data to Characterize Adverse Drug Reaction Quantitatively: Mixed Methods Study.

BACKGROUND: Adverse drug reactions (ADRs), which are the phenotypic manifestations of clinical drug toxicity in humans, are a major concern in precision clinical medicine. A comprehensive evaluation of ADRs is helpful for unbiased supervision of marketed drugs and for discovering new drugs with high success rates. OBJECTIVE: In current practice, drug safety evaluation is often oversimplified to the occurrence or nonoccurrence of ADRs. Given the limitations of current qualitative methods, there is an urgent need for a quantitative evaluation model to improve pharmacovigilance and the accurate assessment of drug safety. METHODS: In this study, we developed a mathematical model, namely the Adverse Drug Reaction Classification System (ADReCS) severity-grading model, for the quantitative characterization of ADR severity, a crucial feature for evaluating the impact of ADRs on human health. The model was constructed by mining millions of real-world historical adverse drug event reports. A new parameter called Severity_score was introduced to measure the severity of ADRs, and upper and lower score boundaries were determined for 5 severity grades. RESULTS: The ADReCS severity-grading model exhibited excellent consistency (99.22%) with the expert-grading system, the Common Terminology Criteria for Adverse Events. Hence, we graded the severity of 6277 standard ADRs for 129,407 drug-ADR pairs. Moreover, we calculated the occurrence rates of 6272 distinct ADRs for 127,763 drug-ADR pairs in large patient populations by mining real-world medication prescriptions. With the quantitative features, we demonstrated example applications in systematically elucidating ADR mechanisms and thereby discovered a list of drugs with improper dosages. CONCLUSIONS: In summary, this study represents the first comprehensive determination of both ADR severity grades and ADR frequencies. This endeavor establishes a strong foundation for future artificial intelligence applications in discovering new drugs with high efficacy and low toxicity. It also heralds a paradigm shift in clinical toxicity research, moving from qualitative description to quantitative evaluation.

Loss-of-function missense variant of AKAP4 induced male infertility through reduced interaction with QRICH2 during sperm flagella development.

Sperm fibrous sheath (FS) is closely related to sperm maturation, capacitation and motility, and A-kinase anchor protein 4 (AKAP4) is the most abundant protein in sperm FS. Previous studies found incomplete sperm FSs and abnormal flagella in Akap4 knockout mice. Meanwhile, it was reported that the partial deletion in AKAP4 is highly relevant to the dysplasia of the FS in an infertile man, and so far, there is no report about male infertility caused by hemizygous AKAP4 variant. Furthermore, the specific mechanisms of how the variant is relevant to the phenotype remain elusive. In this study, we investigated three multiple morphological abnormalities of the sperm flagella-affected men from three independent families (including one consanguine family) carried hemizygous c.C1285T variant in AKAP4. The patients carried this variant, which showed dysplastic sperm FS, and the protein expression of AKAP4 was decreased in flagella, which was further confirmed in HEK-293T cells in vitro. In addition, the co-localization and interaction between AKAP4 and glutamine-rich protein 2 (QRICH2) on the molecular level were identified by immunofluorescence and co-immunoprecipitation (CO-IP). The hemizygous c.1285C > T variant in AKAP4 induced decreased protein expression of QRICH2 in spermatozoa. These results suggested that the normal expression of AKAP4 is required for maintaining the expression of QRICH2 and the decreased protein expression of AKAP4 and QRICH2，as well as the interaction between them induced by the hemizygous variant of AKAP4 caused dysplastic fibrous sheath, which eventually led to reduced sperm motility and male infertility.

PTP-MEG2 regulates quantal size and fusion pore opening through two distinct structural bases and substrates.

Tyrosine phosphorylation of secretion machinery proteins is a crucial regulatory mechanism for exocytosis. However, the participation of protein tyrosine phosphatases (PTPs) in different exocytosis stages has not been defined. Here we demonstrate that PTP-MEG2 controls multiple steps of catecholamine secretion. Biochemical and crystallographic analyses reveal key residues that govern the interaction between PTP-MEG2 and its substrate, a peptide containing the phosphorylated NSF-pY83 site, specify PTP-MEG2 substrate selectivity, and modulate the fusion of catecholamine-containing vesicles. Unexpectedly, delineation of PTP-MEG2 mutants along with the NSF binding interface reveals that PTP-MEG2 controls the fusion pore opening through NSF independent mechanisms. Utilizing bioinformatics search and biochemical and electrochemical screening approaches, we uncover that PTP-MEG2 regulates the opening and extension of the fusion pore by dephosphorylating the DYNAMIN2-pY125 and MUNC18-1-pY145 sites. Further structural and biochemical analyses confirmed the interaction of PTP-MEG2 with MUNC18-1-pY145 or DYNAMIN2-pY125 through a distinct structural basis compared with that of the NSF-pY83 site. Our studies thus provide mechanistic insights in complex exocytosis processes.

The novel BRDT inhibitor NHWD870 shows potential as a male contraceptive in mice.

Small molecule inhibitors of the bromodomain and extraterminal domain (BET) family proteins have emerged as promising options not only for the treatment of multiple cancers but also for disturbing the process of sperm maturation with potential for use as viable contraceptive targets. In this study, we find that the BET family inhibitor NHWD870 and BRDT can bind well in vitro through bioinformatics software prediction and protein binding inhibition experiments. NHWD870 can produce a good contraceptive effect through animal experiments in vivo, and the fertility can be restored to normal after drug withdrawal. Transcriptomics and proteomics results suggest that NHWD870 affects pathways related to spermatogenesis and maturation, further contributing to the male infertility phenotype. Our results show that NHWD870 can induce a complete and reversible contraceptive effect in mice, which is stronger than that of JQ1 and its synthesized derivatives. This study is expected to eventually lead to clinical trials.

ADReCS-Target: target profiles for aiding drug safety research and application.

Delivering safe and effective therapeutic treatment to patients is one of the grand challenges in modern medicine. However, drug safety research has been progressing slowly in recent years, compared to other fields such as biotechnologies and precision medicine, due to the mechanistic complexity of adverse drug reactions (ADRs). To fill up this gap, we develop a new database, the Adverse Drug Reaction Classification System-Target Profile (ADReCS-Target, http://bioinf.xmu.edu.cn/ADReCS-Target), which provides comprehensive information about ADRs caused by drug interaction with protein, gene and genetic variation. In total, ADReCS-Target includes 66,573 pairwise relations, among which 1710 are protein-ADR associations, 2613 are genetic variation-ADR associations, and 63,298 are gene-ADR associations. In a case study of exploring the mechanism of rash, we find that HLAs, C1QA and APOA1 are the key gene players and thus can be potential targets (or biomarkers) in monitoring or countermining rashes. In summary, ADReCS-Target can be a useful resource for the biomedical scientific community by serving researchers in the fields of drug development, clinical pharmacology, precision medicine, and from web lab to high-throughput computational platform. Particularly, it helps to identify drug with better ADR profile and design safer drug therapy regimen.

Biophysical characterization and ligand-binding properties of the elongation factor Tu from Mycobacterium tuberculosis.

Mycobacterium tuberculosis (Mtb) is the key devastating bacterial pathogen responsible for tuberculosis. Increasing emergence of multi-drug-resistant, extensively drug-resistant, and rifampicin/isoniazid-resistant strains of Mtb makes the discovery of validated drug targets an urgent priority. As a vital translational component of the protein biosynthesis system, elongation factor Tu (EF-Tu) is an important molecular switch responsible for selection and binding of the cognate aminoacyl-tRNA to the acceptor site on the ribosome. In addition, EF-Tu from Mtb (MtbEF-Tu) is involved in the initial step of trans-translation which is an effective system for rescuing the stalled ribosomes from non-stop translation complexes under stress conditions. Given its crucial role in protein biosynthesis, EF-Tu is identified as an excellent molecular target for drug design. Here, we reported the recombinant expression, purification, biophysical characterization, and structural modeling of the MtbEF-Tu protein. Our results demonstrated that prokaryotic expression plasmids of pET28a-MtbEF-Tu could be expressed efficiently in Escherichia coli. We successfully purified the 6× His-tagged proteins with a yield of 16.8 mg from 1 l of Luria Bertani medium. Dynamic light scattering experiments showed that MtbEF-Tu existed in a monomeric form, and circular dichroism experiments indicated that MtbEF-Tu was well structured. Moreover, isothermal titration calorimetry experiments displayed that the purified MtbEF-Tu protein possessed intermediate binding affinities for guanosine-5'-triphosphate (GTP) and GDP. The GTP/GDP-binding sites were predicted by flexible molecular docking approach which reveals that GTP/GDP binds to MtbEF-Tu mainly through hydrogen bonds. Our work lays the essential basis for further structural and functional studies of MtbEF-Tu as well as MtbEF-Tu-related novel drug developments.

Evolutionary relationships between seryl-histidine dipeptide and modern serine proteases from the analysis based on mass spectrometry and bioinformatics.

Seryl-histidine dipeptide (Ser-His) has been recognized as the shortest peptide with hydrolysis cleavage activity; however, its protein cleavage spectrum has not yet been fully explored. Here, four differently folded proteins were treated with Ser-His, and the digestion products were evaluated with high-resolution mass spectrometry. The cleavage efficiency and cleavage propensity of Ser-His against these protein substrates were calculated at both the primary and secondary sequence levels. The above experiments show that Ser-His cleaves a broad spectrum of substrate proteins of varying secondary structures. Moreover, Ser-His could cleave at all 20 amino acids with different efficiencies according to the protein, which means that Ser-His has the original digestion function of serine proteases. Furthermore, we collected and compared the catalytic sites and cleavage sites of 340 extant serine proteases derived from 17 representative organisms. A consensus motif Ser-[X]-His was identified as the major pattern at the catalytic sites of serine proteases from all of the organisms represented except Danio rerio, which uses Ser-Lys instead. This finding indicates that Ser-His is the core component of the serine protease catalytic site. Moreover, our analysis revealed that the cleavage sites of modern serine proteases have become more specific over the evolutionary history of this family. Based on the above analysis results, it could be found that Ser-His is likely the original serine protease and maybe the evolutionary core of modern serine proteases.

Competitive Inhibition of Lysine Acetyltransferase 2B by a Small Motif of the Adenoviral Oncoprotein E1A.

The adenovirus early region 1A (E1A) oncoprotein hijacks host cells via direct interactions with many key cellular proteins, such as KAT2B, also known as PCAF (p300/CBP associated factor). E1A binds the histone acetyltransferase (HAT) domain of KAT2B to repress its transcriptional activation. However, the molecular mechanism by which E1A inhibits the HAT activity is not known. Here we demonstrate that a short and relatively conserved N-terminal motif (cNM) in the intrinsically disordered E1A protein is crucial for KAT2B interaction, and inhibits its HAT activity through a direct competition with acetyl-CoA, but not its substrate histone H3. Molecular modeling together with a series of mutagenesis experiments suggests that the major helix of E1A cNM binds to a surface of the acetyl-CoA pocket of the KAT2B HAT domain. Moreover, transient expression of the cNM peptide is sufficient to inhibit KAT2B-specific H3 acetylation H3K14ac in vivo Together, our data define an essential motif cNM in N-terminal E1A as an acetyl-CoA entry blocker that directly associates with the entrance of acetyl-CoA binding pocket to block the HAT domain access to its cofactor.

Analysis of the chemical toxicity effects using the enrichment of Gene Ontology terms and KEGG pathways.

BACKGROUND: Chemical toxicity is one of the major barriers for designing and detecting new chemical entities during drug discovery. Unexpected toxicity of an approved drug may lead to withdrawal from the market and significant loss of the associated costs. Better understanding of the mechanisms underlying various toxicity effects can help eliminate unqualified candidate drugs in early stages, allowing researchers to focus their attention on other more viable candidates. METHODS: In this study, we aimed to understand the mechanisms underlying several toxicity effects using Gene Ontology (GO) terms and KEGG pathways. GO term and KEGG pathway enrichment theories were adopted to encode each chemical, and the minimum redundancy maximum relevance (mRMR) was used to analyze the GO terms and the KEGG pathways. Based on the feature list obtained by the mRMR method, the most related GO terms and KEGG pathways were extracted. RESULTS: Some important GO terms and KEGG pathways were uncovered, which were concluded to be significant for determining chemical toxicity effects. CONCLUSIONS: Several GO terms and KEGG pathways are highly related to all investigated toxicity effects, while some are specific to a certain toxicity effect. GENERAL SIGNIFICANCE: The findings in this study have the potential to further our understanding of different chemical toxicity mechanisms and to assist scientists in developing new chemical toxicity prediction algorithms. This article is part of a Special Issue entitled "System Genetics" Guest Editor: Dr. Yudong Cai and Dr. Tao Huang.

ADReCS: an ontology database for aiding standardization and hierarchical classification of adverse drug reaction terms.

Adverse drug reactions (ADRs) are noxious and unexpected effects during normal drug therapy. They have caused significant clinical burden and been responsible for a large portion of new drug development failure. Molecular understanding and in silico evaluation of drug (or candidate) safety in laboratory is thus so desired, and unfortunately has been largely hindered by misuse of ADR terms. The growing impact of bioinformatics and systems biology in toxicological research also requires a specialized ADR term system that works beyond a simple glossary. Adverse Drug Reaction Classification System (ADReCS; http://bioinf.xmu.edu.cn/ADReCS) is a comprehensive ADR ontology database that provides not only ADR standardization but also hierarchical classification of ADR terms. The ADR terms were pre-assigned with unique digital IDs and at the same time were well organized into a four-level ADR hierarchy tree for building an ADR-ADR relation. Currently, the database covers 6544 standard ADR terms and 34,796 synonyms. It also incorporates information of 1355 single active ingredient drugs and 134,022 drug-ADR pairs. In summary, ADReCS offers an opportunity for direct computation on ADR terms and also provides clues to mining common features underlying ADRs.

Briefing in family characteristics of microRNAs and their applications in cancer research.

MicroRNAs (miRNAs) are endogenous, short, non-coding RNA molecules that are directly involved in the post-transcriptional regulation of gene expression. Dysregulation of miRNAs is usually associated with diseases. Since miRNAs in a family intend to have common functional characteristics, proper assignment of miRNA family becomes heuristic for better understanding of miRNA nature and their potentials in clinic. In this review, we will briefly discuss the recent progress in miRNA research, particularly its impact on protein and its clinical application in cancer research in a view of miRNA family. This article is part of a Special Issue entitled: Computational Proteomics, Systems Biology & Clinical Implications. Guest Editor: Yudong Cai.

High-throughput identification of off-targets for the mechanistic study of severe adverse drug reactions induced by analgesics.

Drugs may induce adverse drug reactions (ADRs) when they unexpectedly bind to proteins other than their therapeutic targets. Identification of these undesired protein binding partners, called off-targets, can facilitate toxicity assessment in the early stages of drug development. In this study, a computational framework was introduced for the exploration of idiosyncratic mechanisms underlying analgesic-induced severe adverse drug reactions (SADRs). The putative analgesic-target interactions were predicted by performing reverse docking of analgesics or their active metabolites against human/mammal protein structures in a high-throughput manner. Subsequently, bioinformatics analyses were undertaken to identify ADR-associated proteins (ADRAPs) and pathways. Using the pathways and ADRAPs that this analysis identified, the mechanisms of SADRs such as cardiac disorders were explored. For instance, 53 putative ADRAPs and 24 pathways were linked with cardiac disorders, of which 10 ADRAPs were confirmed by previous experiments. Moreover, it was inferred that pathways such as base excision repair, glycolysis/glyconeogenesis, ErbB signaling, calcium signaling, and phosphatidyl inositol signaling likely play pivotal roles in drug-induced cardiac disorders. In conclusion, our framework offers an opportunity to globally understand SADRs at the molecular level, which has been difficult to realize through experiments. It also provides some valuable clues for drug repurposing.

PhosMap: An ensemble bioinformatic platform to empower interactive analysis of quantitative phosphoproteomics.

BACKGROUND: Liquid chromatography-mass spectrometry (LC-MS)-based quantitative phosphoproteomics has been widely used to detect thousands of protein phosphorylation modifications simultaneously from the biological specimens. However, the complicated procedures for analyzing phosphoproteomics data has become a bottleneck to widening its application. METHODS: Here, we develop PhosMap, a versatile and scalable tool to accomplish phosphoproteomics data analysis. A standardized phosphorylation data format was created for data analyses, from data preprocessing to downstream bioinformatic analyses such as dimension reduction, differential phosphorylation analysis, kinase activity, survival analysis, and so on. For better usability, we distribute PhosMap as a Docker image for easy local deployment upon any of Windows, Linux, and Mac system. RESULTS: The source code is deposited at https://github.com/BADD-XMU/PhosMap. A free PhosMap webserver (https://huggingface.co/spaces/Bio-Add/PhosMap), with easy-to-follow fashion of dashboards, is curated for interactive data analysis. CONCLUSIONS: PhosMap fills the technical gap of large-scale phosphorylation research by empowering researchers to process their own phosphoproteomics data expediently and efficiently, and facilitates better data interpretation.

PaGeFinder: quantitative identification of spatiotemporal pattern genes.

UNLABELLED: Pattern Gene Finder (PaGeFinder) is a web-based server for on-line detection of gene expression patterns from serial transcriptomic data generated by high-throughput technologies like microarray or next-generation sequencing. Three particular parameters, the specificity measure, the dispersion measure and the contribution measure, were introduced and implemented in PaGeFinder to help quantitative and interactive identification of pattern genes like housekeeping genes, specific (selective) genes and repressed genes. Besides the on-line computation service, the PaGeFinder also provides downloadable Java programs for local detection of gene expression patterns. AVAILABILITY: http://bioinf.xmu.edu.cn:8080/PaGeFinder/index.jsp

The GH18 family of chitinases: their domain architectures, functions and evolutions.

The glycoside hydrolase 18 (GH18) family of chitinases is a multigene family that plays various roles, such as ecdysis, embryonic development, allergic inflammation and so on. Efforts are still needed to reveal their functional diversification in an evolutionary and systematic manner. We collected 85 GH18 genes from eukaryotic representatives. The domain architectures of GH18 proteins were analyzed and several conserved patterns were identified. It was observed that some (11 proteins) GH18 members in Ecdysozoa or fungi possess repeats of catalytic domains and/or chitin-binding domains (ChtBs). The domain repeats are likely to meet requirements for higher efficiency of chitin degradation in chitin-containing species. On the contrary, all vertebrate GH18 proteins contain no more than one catalytic domain or ChtB. The results from homologous analysis, domain architectures, exon arrangements and synteny loci supported two evolutionary paths for the GH18 family. One path experienced gene expansion and contraction several times during evolution, covering most of GH18 members except CHID1 (stabilin-1 interacting partner) and its homologs. Proteins in this path underwent frequent domain gain and loss, as well as domain recombination, that could achieve versatility in function. The other path is comparatively conserved. The CHID1 gene evolved without gene duplication except in Danio rerio. Domain architectures of CHID1 orthologs are all identical. The diverse phylogeny of the GH18 family in arthropod is also presented.