In silico and in vitro study of FLT3 inhibitors and their application in acute myeloid leukemia.

Acute myeloid leukemia (AML) is the most common hematological cancer in the adult population worldwide. Approximately 35% of patients with AML present internal tandem duplication (ITD) mutations in the FMS­like tyrosine kinase 3 (FLT3) receptor associated with poor prognosis, and thus, this receptor is a relevant target for potential therapeutics. Tyrosine kinase inhibitors (TKIs) are used to treat AML; however, their molecular interactions and effects on leukemic cells are poorly understood. The present study aimed to gain insights into the molecular interactions and affinity forces of four TKI drugs (sorafenib, midostaurin, gilteritinib and quizartinib) with the wild­type (WT)­FLT3 and ITD­mutated (ITD­FLT3) structural models of FLT3, in its inactive aspartic acid­phenylalanine­glycine motif (DFG­out) and active aspartic acid­phenylalanine­glycine motif (DFG­in) conformations. Furthermore, the present study evaluated the effects of the second­generation TKIs gilteritinib and quizartinib on cancer cell viability, apoptosis and proliferation in the MV4­11 (ITD­FLT3) and HL60 (WT­FLT3) AML cell lines. Peripheral blood mononuclear cells (PBMCs) from a healthy volunteer were included as an FLT3­negative group. Molecular docking analysis indicated higher affinities of second­generation TKIs for WT­FLT3/DFG­out and WT­FLT3/DFG­in compared with those of the first­generation TKIs. However, the ITD mutation changed the affinity of all TKIs. The in vitro data supported the in silico predictions: MV4­11 cells presented high selective sensibility to gilteritinib and quizartinib compared with the HL60 cells, whereas the drugs had no effect on PBMCs. Thus, the current study presented novel information about molecular interactions between the FLT3 receptors (WT or ITD­mutated) and some of their inhibitors. It also paves the way for the search for novel inhibitory molecules with potential use against AML.

SHARK enables sensitive detection of evolutionary homologs and functional analogs in unalignable and disordered sequences.

Intrinsically disordered regions (IDRs) are structurally flexible protein segments with regulatory functions in multiple contexts, such as in the assembly of biomolecular condensates. Since IDRs undergo more rapid evolution than ordered regions, identifying homology of such poorly conserved regions remains challenging for state-of-the-art alignment-based methods that rely on position-specific conservation of residues. Thus, systematic functional annotation and evolutionary analysis of IDRs have been limited, despite them comprising ~21% of proteins. To accurately assess homology between unalignable sequences, we developed an alignment-free sequence comparison algorithm, SHARK (Similarity/Homology Assessment by Relating K-mers). We trained SHARK-dive, a machine learning homology classifier, which achieved superior performance to standard alignment-based approaches in assessing evolutionary homology in unalignable sequences. Furthermore, it correctly identified dissimilar but functionally analogous IDRs in IDR-replacement experiments reported in the literature, whereas alignment-based tools were incapable of detecting such functional relationships. SHARK-dive not only predicts functionally similar IDRs at a proteome-wide scale but also identifies cryptic sequence properties and motifs that drive remote homology and analogy, thereby providing interpretable and experimentally verifiable hypotheses of the sequence determinants that underlie such relationships. SHARK-dive acts as an alternative to alignment to facilitate systematic analysis and functional annotation of the unalignable protein universe.

Homology-independent targeted insertion-mediated derivation of M1-biased macrophages harbouring Megf10 and CD3ζ from human pluripotent stem cells.

BACKGROUND: Macrophages engineered with chimeric antigen receptors (CAR) are suitable for immunotherapy based on their immunomodulatory activity and ability to infiltrate solid tumours. However, the production and application of genetically edited, highly effective, and mass-produced CAR-modified macrophages (CAR-Ms) are challenging. METHODS: Here, we used homology-independent targeted insertion (HITI) for site-directed CAR integration into the safe-harbour region of human pluripotent stem cells (hPSCs). This approach, together with a simple differentiation protocol, produced stable and highly effective CAR-Ms without heterogeneity. FINDINGS: These engineered cells phagocytosed cancer cells, leading to significant inhibition of cancer-cell proliferation in vitro and in vivo. Furthermore, the engineered CARs, which incorporated a combination of CD3ζ and Megf10 (referred to as FRP5Mζ), markedly enhanced the antitumour effect of CAR-Ms by promoting M1, but not M2, polarisation. FRP5Mζ promoted M1 polarisation via nuclear factor kappa B (NF-κB), ERK, and STAT1 signalling, and concurrently inhibited STAT3 signalling even under M2 conditions. These features of CAR-Ms modulated the tumour microenvironment by activating inflammatory signalling, inducing M1 polarisation of bystander non-CAR macrophages, and enhancing the infiltration of T cells in cancer spheroids. INTERPRETATION: Our findings suggest that CAR-Ms have promise as immunotherapeutics. In conclusion, the guided insertion of CAR containing CD3ζ and Megf10 domains is an effective strategy for the immunotherapy of solid tumours. FUNDING: This work was supported by KRIBB Research Initiative Program Grant (KGM4562431, KGM5282423) and a Korean Fund for Regenerative Medicine (KFRM) grant funded by the Korean government (Ministry of Science and ICT,Ministry of Health and Welfare) (22A0304L1-01).

Precision and efficacy of RNA-guided DNA integration in high-expressing muscle loci.

Gene replacement therapies primarily rely on adeno-associated virus (AAV) vectors for transgene expression. However, episomal expression can decline over time due to vector loss or epigenetic silencing. CRISPR-based integration methods offer promise for long-term transgene insertion. While the development of transgene integration methods has made substantial progress, identifying optimal insertion loci remains challenging. Skeletal muscle is a promising tissue for gene replacement owing to low invasiveness of intramuscular injections, relative proportion of body mass, the multinucleated nature of muscle, and the potential for reduced adverse effects. Leveraging endogenous promoters in skeletal muscle, we evaluated two highly expressing loci using homology-independent targeted integration (HITI) to integrate reporter or therapeutic genes in mouse myoblasts and skeletal muscle tissue. We hijacked the muscle creatine kinase (Ckm) and myoglobin (Mb) promoters by co-delivering CRISPR-Cas9 and a donor plasmid with promoterless constructs encoding green fluorescent protein (GFP) or human Factor IX (hFIX). Additionally, we deeply profiled our genome and transcriptome outcomes from targeted integration and evaluated the safety of the proposed sites. This study introduces a proof-of-concept technology for achieving high-level therapeutic gene expression in skeletal muscle, with potential applications in targeted integration-based medicine and synthetic biology.

Functional Analysis of BARD1 Missense Variants on Homology-Directed Repair in Ovarian and Breast Cancers.

Women with germline BRCA1 mutations face an increased risk of developing breast and ovarian cancers. BARD1 (BRCA1 associated RING domain 1) is an essential heterodimeric partner of BRCA1, and mutations in BARD1 are also associated with these cancers. While BARD1 mutations are recognized for their cancer susceptibility, the exact roles of numerous BARD1 missense mutations remain unclear. In this study, we conducted functional assays to assess the homology-directed DNA repair (HDR) activity of all BARD1 missense substitutions identified in 55 breast and ovarian cancer samples, using the real-world data from the COSMIC and cBioPortal databases. Seven BARD1 variants (V85M, P187A, G491R, R565C, P669L, T719R, and Q730L) were confirmed to impair DNA damage repair. Furthermore, cells harboring these BARD1 variants exhibited increased sensitivity to the chemotherapeutic drugs, cisplatin, and olaparib, compared to cells expressing wild-type BARD1. These findings collectively suggest that these seven missense BARD1 variants are likely pathogenic and may respond well to cisplatin-olaparib combination therapy. This study not only enhances our understanding of BARD1's role in DNA damage repair but also offers valuable insights into predicting therapy responses in patients with specific BARD1 missense mutations.

Unlocking Hope: Paving the Way for a Cutting-Edge Multi-Epitope Dengue Virus Vaccine.

Dengue fever is a significant health issue in Pakistan, demanding a vaccine effective against all the viral strains. This study employs reverse vaccinology to develop potential dengue vaccine candidates (DVAX I-III). The study thoroughly examined conserved areas of dengue virus serotypes 1-4's structural and non-structural proteins. Key viral proteins were analyzed to find antigenic peptides, which were incorporated into vaccine candidates and potentiated with adjuvants. Computational methods predicted peptide structures and evaluated their binding to immune receptors TLR 2, TLR 4, HLA *A1101, and DRB*401. A molecular dynamics simulation lasting 100 ns of the DVAX II-TLR4 complex at different time intervals clearly indicated that the ligand is attached to the receptor. Normal mode analysis assessed the stability and flexibility of these interactions. Encouragingly, all three vaccine candidates demonstrated favorable interactions with these immune receptors and the potential to induce a robust immune response. These findings suggest their safety and warrant further in vivo studies to evaluate their efficacy for clinical development.

Osteocalcin binds to a GPRC6A Venus fly trap allosteric site to positively modulate GPRC6A signaling.

GPRC6A, a member of the Family C G-protein coupled receptors, regulates energy metabolism and sex hormone production and is activated by diverse ligands, including cations, L-amino acids, the osteocalcin (Ocn) peptide and the steroid hormone testosterone. We sought a structural framework for the ability of multiple distinct classes of ligands to active GPRC6A. We created a structural model of GPRC6A using Alphafold2. Using this model we explored a putative orthosteric ligand binding site in the bilobed Venus fly trap (VFT) domain of GPRC6A and two positive allosteric modulator (PAM) sites, one in the VFT and the other in the 7 transmembrane (7TM) domain. We provide evidence that Ocn peptides act as a PAM for GPRC6A by binding to a site in the VFT that is distinct from the orthosteric site for calcium and L-amino acids. In agreement with this prediction, alternatively spliced GPRC6A isoforms 2 and 3, which lack regions of the VFT, and mutations in the computationally predicted Ocn binding site, K352E and H355P, prevent Ocn activation of GPRC6A. These observations explain how dissimilar ligands activate GPRC6A and set the stage to develop novel molecules to activate and inhibit this previously poorly understood receptor.

CHATGPT FOR COMPUTATIONAL TOPOLOGY.

ChatGPT represents a significant milestone in the field of artificial intelligence (AI), finding widespread applications across diverse domains. However, its effectiveness in mathematical contexts has been somewhat constrained by its susceptibility to conceptual errors. Concurrently, topological data analysis (TDA), a relatively new discipline, has garnered substantial interest in recent years. Nonetheless, the advancement of TDA is impeded by the limited understanding of computational algorithms and coding proficiency among theoreticians. This work endeavors to bridge the gap between theoretical topological concepts and their practical implementation in computational topology through the utilization of ChatGPT. We showcase how a pure theoretician, devoid of computational experience and coding skills, can effectively transform mathematical formulations and concepts into functional codes for computational topology with the assistance of ChatGPT. Our strategy outlines a productive process wherein a mathematician trains ChatGPT on pure mathematical concepts, steers ChatGPT towards generating computational topology codes, and subsequently validates the generated codes using established examples. Our specific case studies encompass the computation of Betti numbers, Laplacian matrices, and Dirac matrices for simplicial complexes, as well as the persistence of various homologies and Laplacians. Furthermore, we explore the application of ChatGPT in computing recently developed topological theories for hypergraphs and digraphs, as well as the persistent harmonic space, which has not been computed in the literature, to the best of our knowledge. This work serves as an initial step towards effectively transforming pure mathematical theories into practical computational tools, with the ultimate goal of enabling real applications across diverse fields.

Adapting to change: resolving the dynamic and dual roles of NCK1 and NCK2.

Adaptor proteins play central roles in the assembly of molecular complexes and co-ordinated activation of specific pathways. Through their modular domain structure, the NCK family of adaptor proteins (NCK1 and NCK2) link protein targets via their single SRC Homology (SH) 2 and three SH3 domains. Classically, their SH2 domain binds to phosphotyrosine motif-containing receptors (e.g. receptor tyrosine kinases), while their SH3 domains bind polyproline motif-containing cytoplasmic effectors. Due to these functions being established for both NCK1 and NCK2, their roles were inaccurately assumed to be redundant. However, in contrast with this previously held view, NCK1 and NCK2 now have a growing list of paralog-specific functions, which underscores the need to further explore their differences. Here we review current evidence detailing how these two paralogs are unique, including differences in their gene/protein regulation, binding partners and overall contributions to cellular functions. To help explain these contrasting characteristics, we then discuss SH2/SH3 structural features, disordered interdomain linker regions and post-translational modifications. Together, this review seeks to highlight the importance of distinguishing NCK1 and NCK2 in research and to pave the way for investigations into the origins of their interaction specificity.

A stomate by any other name? The open question of hornwort gametophytic pores, their homology, and implications for the evolution of stomates.

Advances in bryophyte genomics and the phylogenetic recovery of hornworts, mosses, and liverworts as a clade have spurred considerable recent interest in character evolution among early embryophytes. Discussion of stomatal evolution, however, has been incomplete; the result of the neglect of certain potential stomate homologues, namely the two-celled epidermal gametophytic pores of hornworts (typically referred to as 'mucilage clefts'). Confusion over the potential homology of these structures is the consequence of a relatively recent consensus that hornwort gametophytic pores ('HGPs' - our term) are not homologous to stomates. We explore the occurrence and diverse functions of stomates throughout the evolutionary history and diversity of extinct and extant embryophytes. We then address arguments for and against homology between known sporophyte- and gametophyte-borne stomates and HGPs and conclude that there is little to no evidence that contradicts the hypothesis of homology. We propose that 'intergenerational heterotopy' might well account for the novel expression of stomates in gametophytes of hornworts, if stomates first evolved in the sporophyte generation of embryophytes. We then explore phylogenetically based hypotheses for the evolution of stomates in both the gametophyte and sporophyte generations of early lineages of embryophytes.

PungentDB: Bridging traditional Chinese medicine of medicine food homology and modern food flavor chemistry.

Merging traditional Chinese medicine's (TCM) principles of medicine-food homology with modern flavor chemistry, this research unveils PungentDB (http://www.pungentdb.org.cn/home), a database documenting 205 unique pungent flavor compounds from 231 TCMs. It provides detailed insights into their chemical attributes, biological targets (including IC50/EC50 values), and molecular structures (2D/3D), enriched with visualizations of target organ distribution and protein structures, exploring the pungent flavor space with the help of a feature-rich visual interface. This collection, derived from over 3249 sources and highlighting 9129 targets, delves into the compounds' unique pungent flavors-taste, aroma, and thermal sensations-and their interaction with taste and olfactory receptors. PungentDB bridges ancient wisdom and culinary innovation, offering a nuanced exploration of pungent flavors' role in enhancing food quality, safety, and sensory experiences. This initiative propels flavor chemistry forward, serving as a pivotal resource for food science advancement and the innovative application of pungent flavors.

Discovery of a novel mutation F184S (c.551T>C) in GATA4 gene causing congenital heart disease in a consanguineous Saudi family.

Background & aim: Congenital heart disease (CHD) is the most common cause of non-infectious deaths in infants worldwide. However, the molecular mechanisms underlying CHD remain unclear. Approximately 30 % of the causes are believed to be genetic mutations and chromosomal abnormalities. In this study, we aimed to identify the genetic causes of CHD in consanguineous families. Methods: Fourth-generation pedigrees with CHD were recruited. The main cardiac features of the patient included absence of the right pulmonary artery and a large dilated left pulmonary artery. To determine the underlying genetic cause, whole-exome sequencing was performed and subsequently confirmed using Sanger sequencing and different online databases to study the pathogenesis of the identified gene mutation. An in-silico homology model was created using the Alphafold homology model structure of GATA4 (AF-P43694-F1). The missense3D online program was used to evaluate the structural alterations. Results: We identified a deleterious mutation c.551T > C (p.Phe184Ser) in GATA4. GATA4 is a highly conserved zinc-finger transcription factor, and its continuous expression is essential for cardiogenesis during embryogenesis. The in-silico model suggested a compromised binding efficiency with other proteins. Several variant interpretation algorithms indicated that the F184S missense variant in GATA4 is damaging, whereas HOPE analysis indicated the functional impairment of DNA binding of transcription factors and zinc-ion binding activities of GATA4. Conclusion: The variant identified in GATA4 appears to cause recessive CHD in the family. In silico analysis suggested that this variant was damaging and caused multiple structural and functional aberrations. This study may support prenatal screening of the fetus in this family to prevent diseases in new generations.

Adiponectin Resistance in Obesity: Adiponectin Leptin/Insulin Interaction.

The adiponectin (APN) levels in obesity are negatively correlated with chronic subclinical inflammation markers. The hypertrophic adipocytes cause obesity-linked insulin resistance and metabolic syndrome. Furthermore, macrophage polarization is a key determinant regulating adiponectin receptor (AdipoR1/R2) expression and differential adiponectin-mediated macrophage inflammatory responses in obese individuals. In addition to decrease in adiponectin concentrations, the decline in AdipoR1/R2 messenger ribonucleic acid (mRNA) expression leads to a decrement in adiponectin binding to cell membrane, and this turns into attenuation in the adiponectin effects. This is defined as APN resistance, and it is linked with insulin resistance in high-fat diet-fed subjects. The insulin-resistant group has a significantly higher leptin-to-APN ratio. The leptin-to-APN ratio is more than twofold higher in obese individuals. An increase in expression of AdipoRs restores insulin sensitivity and ß-oxidation of fatty acids via triggering intracellular signal cascades. The ratio of high molecular weight to total APN is defined as the APN sensitivity index (ASI). This index is correlated to insulin sensitivity. Homeostasis model of assessment (HOMA)-APN and HOMA-estimated insulin resistance (HOMA-IR) are the most suitable methods to estimate the metabolic risk in metabolic syndrome. While morbidly obese patients display a significantly higher plasma leptin and soluble (s)E-selectin concentrations, leptin-to-APN ratio, there is a significant negative correlation between leptin-to-APN ratio and sP-selectin in obese patients. When comparing the metabolic dysregulated obese group with the metabolically healthy obese group, postprandial triglyceride clearance, insulin resistance, and leptin resistance are significantly delayed following the oral fat tolerance test in the first group. A neuropeptide, Spexin (SPX), is positively correlated with the quantitative insulin sensitivity check index (QUICKI) and APN. APN resistance together with insulin resistance forms a vicious cycle. Despite normal or high APN levels, an impaired post-receptor signaling due to adaptor protein-containing pleckstrin homology domain, phosphotyrosine-binding domain, and leucine zipper motif 1 (APPL1)/APPL2 may alter APN efficiency and activity. However, APPL2 blocks adiponectin signaling through AdipoR1 and AdipoR2 because of the competitive inhibition of APPL1. APPL1, the intracellular binding partner of AdipoRs, is also an important mediator of adiponectin-dependent insulin sensitization. The elevated adiponectin levels with adiponectin resistance are compensatory responses in the condition of an unusual discordance between insulin resistance and APN unresponsiveness. Hypothalamic recombinant adeno-associated virus (rAAV)-leptin (Lep) gene therapy reduces serum APN levels, and it is a more efficient strategy for long-term weight maintenance.

Lonicerae Japonicae Flos with the homology of medicine and food: a review of active ingredients, anticancer mechanisms, pharmacokinetics, quality control, toxicity and applications.

Lonicerae Japonicae Flos (LJF, called Jinyinhua in China), comes from the dried flower buds or flowers to be opened of Lonicera japonica Thunb. in the Lonicera family. It has a long history of medicinal use and has a wide range of application prospects. As modern research advances, an increasing number of scientific experiments have demonstrated the anticancer potential of LJF. However, there is a notable absence of systematic reports detailing the anti-tumor effects of LJF. This review integrates the principles of Traditional Chinese Medicine (TCM) with contemporary pharmacological techniques, drawing upon literature from authoritative databases such as PubMed, CNKI, and WanFang to conduct a comprehensive study of LJF. Notably, a total of 507 compounds have been isolated and characterized from the plant to date, which include volatile oils, organic acids, flavonoids, iridoids, triterpenes and triterpenoid saponins. Pharmacological studies have demonstrated that LJF extract, along with components such as chlorogenic acid, luteolin, rutin, luteoloside, hyperoside and isochlorogenic acid, exhibits potential anticancer activities. Consequently, we have conducted a comprehensive review and summary of the mechanisms of action and clinical applications of these components. Furthermore, we have detailed the pharmacokinetics, quality control, and toxicity of LJF, while also discussing its prospective applications in the fields of biomedicine and preventive healthcare. It is hoped that these studies will provide valuable reference for the clinical research, development, and application of LJF.

Adaptor protein Src-homology 2 domain containing E (SH2E) deficiency induces heart defect in zebrafish.

Adaptor proteins play crucial roles in signal transduction across diverse signaling pathways. Src-homology 2 domain-containing E (SH2E) is the adaptor protein highly expressed in vascular endothelial cells and myocardium during zebrafish embryogenesis. In this study we investigated the function and mechanisms of SH2E in cardiogenesis. We first analyzed the spatiotemporal expression of SH2E and then constructed zebrafish lines with SH2E deficiency using the CRISPR-Cas9 system. We showed that homozygous mutants developed progressive pericardial edema (PCE), dilated atrium, abnormal atrioventricular looping and thickened atrioventricular wall from 3 days post fertilization (dpf) until death; inducible overexpression of SH2E was able to partially rescue the PCE phenotype. Using transcriptome sequencing analysis, we demonstrated that the MAPK/ERK and NF-κB signaling pathways might be involved in SH2E-deficiency-caused PCE. This study underscores the pivotal role of SH2E in cardiogenesis, and might help to identify innovative diagnostic techniques and therapeutic strategies for congenital heart disease.

Inhibitor binding and disruption of coupled motions in MmpL3 protein: Unraveling the mechanism of trehalose monomycolate transport.

Mycobacterial membrane protein Large 3 (MmpL3) of Mycobacterium tuberculosis (Mtb) is crucial for the translocation of trehalose monomycolate (TMM) across the inner bacterial cell membrane, making it a promising target for anti-tuberculosis (TB) drug development. While several structural, microbiological, and in vitro studies have provided significant insights, the precise mechanisms underlying TMM transport by MmpL3 and its inhibition remain incompletely understood at the atomic level. In this study, molecular dynamic (MD) simulations for the apo form and seven inhibitor-bound forms of Mtb MmpL3 were carried out to obtain a thorough comprehension of the protein's dynamics and function. MD simulations revealed that the seven inhibitors in this work stably bind to the central channel of the transmembrane domain and primarily forming hydrogen bonds with ASP251, ASP640, or both residues. Through dynamical cross-correlation matrix and principal component analysis analyses, several types of coupled motions between different domains were observed in the apo state, and distinct conformational states were identified using Markov state model analysis. These coupled motions and varied conformational states likely contribute to the transport of TMM. However, simulations of inhibitor-bound MmpL3 showed an enlargement of the proton channel, potentially disrupting coupled motions. This indicates that inhibitors may impair MmpL3's transport function by directly blocking the proton channel, thereby hindering coordinated domain movements and indirectly affecting TMM translocation.

A novel three-part pharynx and its parallel evolution within symbiotic marine nematodes (Desmodoroidea, Stilbonematinae).

Stilbonematinae are nematodes commonly found in shallow marine sands. They are overgrown by a genus- and species-specific coat of chemoautotrophic sulphur-oxidizing ectosymbiotic bacteria which profit from the vertical migration of their hosts through the chemocline by alternately gaining access to oxidizing and reducing chemical species, while in return, the host feeds on its symbionts. The subfamily exhibits a large morphological variability; e.g. the anterior pharynx is cylindrical in genera possessing a voluminous coat, but species with a bacterial monolayer possess a distinctly swollen corpus and therefore a tripartite pharynx. Since 18S-based phylogenetic analyses do not show close relationships between corpus-bearing species, we investigated the pharynx morphology using phalloidin staining in combination with confocal laser scanning microscopy, transmission electron microscopy and light microscopy in order to assess an independent evolution. The class-wide stable position of the subventral pharynx ampullae was used as a morphological marker. Ampullae are positioned at the anterior-most end of the isthmus in Cyathorobbea, further posterior in Catanema and Robbea and inside the corpus in Laxus oneistus. We therefore conclude an independent evolution of corpus enlargements within Stilbonematinae. This further suggests that pharynx morphology is driven by the volume of the symbiotic bacterial coat rather than phylogeny. Based on an existing mathematical model, an enlarged corpus should enable its bearer to ingest food in smaller quantities, in gourmet style, whereas a cylindrical pharynx would restrict its bearer to ancestral gourmand feeding. A review of pharynx types of Nematoda showed that the Stilbonematinae pharynx is substantially different compared to other tripartite pharynges. The lack of pharyngeal tubes and valves, the undivided corpus and evenly distributed nuclei in the isthmus warrant the definition of the "stilbonematoid" three-part pharynx.

Phylogenetic analyses and antigenic characterization of foot-and-mouth disease virus PanAsia lineage circulating in China between 1999 and 2023.

Foot-and-mouth disease (FMD) is one of the most important transboundary animal diseases caused by foot-and-mouth disease virus (FMDV), leading to significant economic losses worldwide. The first report of PanAsia lineage of FMDV in China was in 1999. Since 2011, 18 outbreaks attributed to PanAsia lineage viruses have been reported across 7 provinces or municipality in China. Phylogenetic analysis indicated that these PanAsia strains were clustered into three distinct clades (clade 1, clade 2, and clade 3), with nucleotide homology ranging from 91.4% to 100%. The outbreaks of FMD caused by clade 1 strains occurred around 1999 when this lineage was prevalent globally. Clade 2 strains dominated from 2011 to 2013, while clade 3 strains were prevalent during 2018-2019, sharing only 93% homology with clade 2 strains and 91% with clade 1 strains. Tracing analysis showed that these outbreaks represented 3 distinct introductions of PanAsia viruses into China. Virus neutralization tests (VNT) have demonstrated that current commercial vaccines are effective to protect susceptible animals against these strains (r1 > 0.3). However, the growing demand for livestock has promoted animal movement and encouraged the exchange of products, services, and materials between countries, thereby heightening the risk of exotic strain incursions. Therefore, it is imperative to reinforce border controls and limit animal movements among various Asian countries continually to reduce the risk of new transboundary diseases, such as FMD incursion. Additionally, PanAsia-2 strains need to be taken seriously to prevent its incursions, and the relevant vaccines against PanAsia-2 strains needs to be stockpiled in preparation for any possible incursion.

Relational Persistent Homology for Multispecies Data with Application to the Tumor Microenvironment.

Topological data analysis (TDA) is an active field of mathematics for quantifying shape in complex data. Standard methods in TDA such as persistent homology (PH) are typically focused on the analysis of data consisting of a single entity (e.g., cells or molecular species). However, state-of-the-art data collection techniques now generate exquisitely detailed multispecies data, prompting a need for methods that can examine and quantify the relations among them. Such heterogeneous data types arise in many contexts, ranging from biomedical imaging, geospatial analysis, to species ecology. Here, we propose two methods for encoding spatial relations among different data types that are based on Dowker complexes and Witness complexes. We apply the methods to synthetic multispecies data of a tumor microenvironment and analyze topological features that capture relations between different cell types, e.g., blood vessels, macrophages, tumor cells, and necrotic cells. We demonstrate that relational topological features can extract biological insight, including the dominant immune cell phenotype (an important predictor of patient prognosis) and the parameter regimes of a data-generating model. The methods provide a quantitative perspective on the relational analysis of multispecies spatial data, overcome the limits of traditional PH, and are readily computable.