Pathobiological signatures of dysbiotic lung injury in pediatric patients undergoing stem cell transplantation.

Hematopoietic cell transplantation (HCT) uses cytotoxic chemotherapy and/or radiation followed by intravenous infusion of stem cells to cure malignancies, bone marrow failure and inborn errors of immunity, hemoglobin and metabolism. Lung injury is a known complication of the process, due in part to disruption in the pulmonary microenvironment by insults such as infection, alloreactive inflammation and cellular toxicity. How microorganisms, immunity and the respiratory epithelium interact to contribute to lung injury is uncertain, limiting the development of prevention and treatment strategies. Here we used 278 bronchoalveolar lavage (BAL) fluid samples to study the lung microenvironment in 229 pediatric patients who have undergone HCT treated at 32 children's hospitals between 2014 and 2022. By leveraging paired microbiome and human gene expression data, we identified high-risk BAL compositions associated with in-hospital mortality (P = 0.007). Disadvantageous profiles included bacterial overgrowth with neutrophilic inflammation, microbiome contraction with epithelial fibroproliferation and profound commensal depletion with viral and staphylococcal enrichment, lymphocytic activation and cellular injury, and were replicated in an independent cohort from the Netherlands (P = 0.022). In addition, a broad array of previously occult pathogens was identified, as well as a strong link between antibiotic exposure, commensal bacterial depletion and enrichment of viruses and fungi. Together these lung-immune system-microorganism interactions clarify the important drivers of fatal lung injury in pediatric patients who have undergone HCT. Further investigation is needed to determine how personalized interpretation of heterogeneous pulmonary microenvironments may be used to improve pediatric HCT outcomes.

Emulation and evaluation of tumor cell combined chemotherapy in isotropic/anisotropic collagen fiber microenvironments.

The rapid emergence of anisotropic collagen fibers in the tissue microenvironment is a critical transition point in late-stage breast cancer. Specifically, the fiber orientation facilitates the likelihood of high-speed tumor cell invasion and metastasis, which pose lethal threats to patients. Thus, based on this transition point, one key issue is how to determine and evaluate efficient combination chemotherapy treatments in late-stage cancer. In this study, we designed a collagen microarray chip containing 241 high-throughput microchambers with embedded metastatic breast cancer cell MDA-MB-231-RFP. By utilizing collagen's unique structure and hydromechanical properties, the chip constructed three-dimensional isotropic and anisotropic collagen fiber structures to emulate the tumor cell microenvironment at early and late stages. We injected different chemotherapeutic drugs into its four channels and obtained composite biochemical concentration profiles. Our results demonstrate that anisotropic collagen fibers promote cell proliferation and migration more than isotropic collagen fibers, suggesting that the geometric arrangement of fibers plays an important role in regulating cell behavior. Moreover, the presence of anisotropic collagen fibers may be a potential factor leading to the poor efficacy of combined chemotherapy in late-stage breast cancer. We investigated the efficacy of various chemotherapy drugs using cell proliferation inhibitors paclitaxel and gemcitabine and tumor cell migration inhibitors 7rh and PP2. To ensure the validity of our findings, we followed a systematic approach that involved testing the inhibitory effects of these drugs. According to our results, the drug combinations' effectiveness could be ordered as follows: paclitaxel + gemcitabine > gemcitabine + 7rh > PP2 + paclitaxel > 7rh + PP2. This study shows that the biomimetic chip system not only facilitates the creation of a realistic in vitro model for examining the cell migration mechanism in late-stage breast cancer but also has the potential to function as an effective tool for future chemotherapy assessment and personalized medicine.

Outcomes from hematopoietic stem cell transplantation following treosulfan-based conditioning: A clinical and pharmacokinetic analysis.

BACKGROUND: The aims of this study are to report our experience with treosulfan-based conditioning regimens for patients with non-malignant hematologic conditions, correlating clinical outcomes at different time points post-transplant with treosulfan exposure (AUC). METHODS: This study was a single-center observational study investigating overall survival (OS), disease-free survival (DFS), and event-free survival (EFS) end-points post-transplant. The consequences of treosulfan AUC with respect to toxicity, correction of underlying disease, and long-term chimerism were also explored using pharmacokinetic analysis. RESULTS: Forty-six patients received 49 transplants with treosulfan and fludarabine-based conditioning between 2005 and 2023. Twenty-four patients also received thiotepa. Donor chimerism was assessed on either whole blood or sorted cell lines at different time points post-transplant. Thirty-nine patients received treosulfan pharmacokinetic assessment to evaluate cumulative AUC, with five infants receiving real-time assessment to facilitate daily dose adjustment. OS, DFS, and EFS were 87%, 81%, and 69%, respectively. Median follow-up was 32.1 months (range 0.82-160 months) following transplant. Lower EFS was associated with patient age (<1 year; p = .057) and lower cumulative treosulfan dose (<42 g/m2; p = .003). Stable donor chimerism in B-cell, NK-cell, and granulocyte lineages at 1-year post-transplant were more prevalent in patients receiving thiotepa conditioning. Two infants required daily dose adjustment to treosulfan to avoid high AUC. CONCLUSIONS: Excellent clinical outcomes and stable chimerism were observed in this patient series. The addition of thiotepa conferred no significant toxicity and trended toward sustained ongoing donor engraftment. Correlating treosulfan AUC with long-term patient outcomes is required.

Langerhans cell histiocytosis: NACHO update on progress, chaos, and opportunity on the path to rational cures.

Langerhans cell histiocytosis (LCH) is a myeloid neoplastic disorder characterized by lesions with CD1a-positive/Langerin (CD207)-positive histiocytes and inflammatory infiltrate that can cause local tissue damage and systemic inflammation. Clinical presentations range from single lesions with minimal impact to life-threatening disseminated disease. Therapy for systemic LCH has been established through serial trials empirically testing different chemotherapy agents and durations of therapy. However, fewer than 50% of patients who have disseminated disease are cured with the current standard-of-care vinblastine/prednisone/(mercaptopurine), and treatment failure is associated with long-term morbidity, including the risk of LCH-associated neurodegeneration. Historically, the nature of LCH-whether a reactive condition versus a neoplastic/malignant condition-was uncertain. Over the past 15 years, seminal discoveries have broadly defined LCH pathogenesis; specifically, activating mitogen-activated protein kinase pathway mutations (most frequently, BRAFV600E) in myeloid precursors drive lesion formation. LCH therefore is a clonal neoplastic disorder, although secondary inflammatory features contribute to the disease. These paradigm-changing insights offer a promise of rational cures for patients based on individual mutations, clonal reservoirs, and extent of disease. However, the pace of clinical trial development behind lags the kinetics of translational discovery. In this review, the authors discuss the current understanding of LCH biology, clinical characteristics, therapeutic strategies, and opportunities to improve outcomes for every patient through coordinated agent prioritization and clinical trial efforts.

Different impact of short-term and long-term hindlimb disuse on bone homeostasis.

Disuse osteoporosis is one of the major problems of bone health which commonly occurs in astronauts during long-term spaceflight and bedridden patients. However, the mechanisms underlying such mechanical unloading induced bone loss have not been fully understood. In this study, we employed hindlimb-unloading mice models with different length of tail suspension to investigate if the bone loss was regulated by distinct factors under different duration of disuse. Our micro-CT results showed more significant decrease of bone mass in 6W (6-week) tail-suspension mice compared to the 1W (1-week) tail-suspension ones, as indicated by greater reduction of BV/TV, Tb.N, B.Ar/T.Ar and Ct.Th. RNA-sequencing results showed significant effects of hindlimb disuse on cell locomotion and immune system process which could cause bone loss.Real-time quantitative PCR results indicated a greater number of bone formation related genes that were downregulated in short-term tail-suspension mice compared to the long-term ones. It is, thus, suggested while sustained hindlimb unloading continuously contributes to bone loss, molecular regulation of bone homeostasis tends to reach a balance during this process.

Fine-mapping analysis including over 254,000 East Asian and European descendants identifies 136 putative colorectal cancer susceptibility genes.

Genome-wide association studies (GWAS) have identified more than 200 common genetic variants independently associated with colorectal cancer (CRC) risk, but the causal variants and target genes are mostly unknown. We sought to fine-map all known CRC risk loci using GWAS data from 100,204 cases and 154,587 controls of East Asian and European ancestry. Our stepwise conditional analyses revealed 238 independent association signals of CRC risk, each with a set of credible causal variants (CCVs), of which 28 signals had a single CCV. Our cis-eQTL/mQTL and colocalization analyses using colorectal tissue-specific transcriptome and methylome data separately from 1299 and 321 individuals, along with functional genomic investigation, uncovered 136 putative CRC susceptibility genes, including 56 genes not previously reported. Analyses of single-cell RNA-seq data from colorectal tissues revealed 17 putative CRC susceptibility genes with distinct expression patterns in specific cell types. Analyses of whole exome sequencing data provided additional support for several target genes identified in this study as CRC susceptibility genes. Enrichment analyses of the 136 genes uncover pathways not previously linked to CRC risk. Our study substantially expanded association signals for CRC and provided additional insight into the biological mechanisms underlying CRC development.

Circulating 25-hydroxyvitamin D and survival outcomes of colorectal cancer: evidence from population-based prospective cohorts and Mendelian randomisation.

BACKGROUND: To investigate the association between circulating 25-hydroxyvitamin D (25-OHD) and colorectal cancer (CRC) survival outcomes. METHODS: We conducted analyses among the Study of Colorectal Cancer in Scotland (SOCCS) and the UK Biobank (UKBB). Both cancer-specific survival (CSS) and overall survival (OS) outcomes were examined. The 25-OHD levels were categorised into three groups, and multi-variable Cox-proportional hazard models were applied to estimate hazard ratios (HRs). We performed individual-level Mendelian randomisation (MR) through the generated polygenic risk scores (PRS) of 25-OHD and summary-level MR using the inverse-variance weighted (IVW) method. RESULTS: We observed significantly poorer CSS (HR = 0.65,95%CI = 0.55-0.76,P = 1.03 × 10-7) and OS (HR = 0.66,95%CI = 0.58-0.75,P = 8.15 × 10-11) in patients with the lowest compared to those with the highest 25-OHD after adjusting for covariates. These associations remained across patients with varied tumour sites and stages. However, we found no significant association between 25-OHD PRS and either CSS (HR = 0.98,95%CI = 0.80-1.19,P = 0.83) or OS (HR = 1.07,95%CI = 0.91-1.25,P = 0.42). Furthermore, we found no evidence for causal effects by conducting summary-level MR analysis for either CSS (IVW:HR = 1.04,95%CI = 0.85-1.28,P = 0.70) or OS (IVW:HR = 1.10,95%CI = 0.93-1.31,P = 0.25). CONCLUSION: This study supports the observed association between lower circulating 25-OHD and poorer survival outcomes for CRC patients. Whilst the genotype-specific association between better outcomes and higher 25-OHD is intriguing, we found no support for causality using MR approaches.

Cardiometabolic biomarker patterns associated with cardiac MRI defined fibrosis and microvascular dysfunction in patients with heart failure with preserved ejection fraction.

Introduction: Heart failure with preserved ejection fraction (HFpEF) is a complex disease process influenced by metabolic disorders, systemic inflammation, myocardial fibrosis, and microvascular dysfunction. The goal of our study is to identify potential relationships between plasma biomarkers and cardiac magnetic resonance (CMR) imaging markers in patients with HFpEF. Methods: Nineteen subjects with HFpEF and 15 age-matched healthy controls were enrolled and underwent multiparametric CMR and plasma biomarker analysis using the Olink® Cardiometabolic Panel (Olink Proteomics, Uppsala, Sweden). Partial least squares discriminant analysis (PLS-DA) was used to characterize CMR and biomarker variables that differentiate the subject groups into two principal components. Orthogonal projection to latent structures by partial least squares (OPLS) analysis was used to identify biomarker patterns that correlate with myocardial perfusion reserve (MPR) and extracellular volume (ECV) mapping. Results: A PLS-DA could differentiate between HFpEF and normal controls with two significant components explaining 79% (Q2 = 0.47) of the differences. For OPLS, there were 7 biomarkers that significantly correlated with ECV (R2 = 0.85, Q = 0.53) and 6 biomarkers that significantly correlated with MPR (R2 = 0.92, Q2 = 0.32). Only 1 biomarker significantly correlated with both ECV and MPR. Discussion: Patients with HFpEF have unique imaging and biomarker patterns that suggest mechanisms associated with metabolic disease, inflammation, fibrosis and microvascular dysfunction.

Peeling off the layers from microbial dark matter (MDM): recent advances, future challenges, and opportunities.

Microbes represent the most common organisms on Earth; however, less than 2% of microbial species in the environment can undergo cultivation for study under laboratory conditions, and the rest of the enigmatic, microbial world remains mysterious, constituting a kind of "microbial dark matter" (MDM). In the last two decades, remarkable progress has been made in culture-dependent and culture-independent techniques. More recently, studies of MDM have relied on culture-independent techniques to recover genetic material through either unicellular genomics or shotgun metagenomics to construct single-amplified genomes (SAGs) and metagenome-assembled genomes (MAGs), respectively, which provide information about evolution and metabolism. Despite the remarkable progress made in the past decades, the functional diversity of MDM still remains uncharacterized. This review comprehensively summarizes the recently developed culture-dependent and culture-independent techniques for characterizing MDM, discussing major challenges, opportunities, and potential applications. These activities contribute to expanding our knowledge of the microbial world and have implications for various fields including Biotechnology, Bioprospecting, Functional genomics, Medicine, Evolutionary and Planetary biology. Overall, this review aims to peel off the layers from MDM, shed light on recent advancements, identify future challenges, and illuminate the exciting opportunities that lie ahead in unraveling the secrets of this intriguing microbial realm.

Exploiting Multimetallic Cooperativity in the Ring-Opening Polymerization of Cyclic Esters and Ethers.

The use of multimetallic complexes is a rapidly advancing route to enhance catalyst performance in the ring-opening polymerization of cyclic esters and ethers. Multimetallic catalysts often outperform their monometallic analogues in terms of reactivity and/or polymerization control, and these improvements are typically attributed to "multimetallic cooperativity". Yet the origins of multimetallic cooperativity often remain unclear. This review explores the key factors underpinning multimetallic cooperativity, including metal-metal distances, the flexibility, electronics and conformation of the ligand framework, and the coordination environment of the metal centers. Emerging trends are discussed to provide insights into why cooperativity occurs and how to harness cooperativity for the development of highly efficient multimetallic catalysts.

Upper and/or Lower Respiratory Tract Infection Caused by Human Metapneumovirus After Allogeneic Hematopoietic Stem Cell Transplantation.

BACKGROUND: Human metapneumovirus (hMPV) epidemiology, clinical characteristics and risk factors for poor outcome after allogeneic stem cell transplantation (allo-HCT) remain a poorly investigated area. METHODS: This retrospective multicenter cohort study examined the epidemiology, clinical characteristics, and risk factors for poor outcomes associated with human metapneumovirus (hMPV) infections in recipients of allo-HCT. RESULTS: We included 428 allo-HCT recipients who developed 438 hMPV infection episodes between January 2012 and January 2019. Most recipients were adults (93%). hMPV infections were diagnosed at a median of 373 days after allo-HCT. The infections were categorized as upper respiratory tract disease (URTD) or lower respiratory tract disease (LRTD), with 60% and 40% of cases, respectively. Patients with hMPV LRTD experienced the infection earlier in the transplant course and had higher rates of lymphopenia, neutropenia, corticosteroid use, and ribavirin therapy. Multivariate analysis identified lymphopenia and corticosteroid use (>30 mg/d) as independent risk factors for LRTD occurrence. The overall mortality at day 30 after hMPV detection was 2% for URTD, 12% for possible LRTD, and 21% for proven LRTD. Lymphopenia was the only independent risk factor associated with day 30 mortality in LRTD cases. CONCLUSIONS: These findings highlight the significance of lymphopenia and corticosteroid use in the development and severity of hMPV infections after allo-HCT, with lymphopenia being a predictor of higher mortality in LRTD cases.

Oxidative stress and antioxidants in cataract development.

PURPOSE OF REVIEW: Oxidative stress plays a central role in cataract pathogenesis, a leading cause of global blindness. This review delves into the role of oxidative stress in cataract development and key biomarkers - glutathione (GSH), superoxide dismutase (SOD), malondialdehyde (MDA), and 4-hydroxynonenal (4-HNE) - to clarify their functions and potential applications in predictive diagnostics and therapies. RECENT FINDINGS: Antioxidants serve as pivotal markers in cataract pathogenesis. GSH affects the central lens due to factors such as enzyme depletion and altered connexin expression, impairing GSH diffusion. Age-related oxidative stress may hinder GSH transport via connexin channels or an internal microcirculation system. N-acetylcysteine, a GSH precursor, shows promise in mitigating lens opacity when applied topically. Additionally, SOD, particularly SOD1, correlates with increased cataract development and gel formulations have exhibited protective effects against posterior subscapular cataracts. Lastly, markers of lipid peroxidation, MDA and 4-HNE, have been shown to reflect disease severity. Studies suggest a potential link between 4-HNE and connexin channel modification, possibly contributing to reduced GSH levels. SUMMARY: Oxidative stress is a significant contributor to cataract development, underscoring the importance of antioxidants in diagnosis and treatment. Notably, GSH depletion, SOD decline, and lipid peroxidation markers are pivotal factors in cataract pathogenesis, offering promising avenues for both diagnosis and therapeutic intervention.

Pulmonary microbiome and transcriptome signatures reveal distinct pathobiologic states associated with mortality in two cohorts of pediatric stem cell transplant patients.

Lung injury is a major determinant of survival after pediatric hematopoietic cell transplantation (HCT). A deeper understanding of the relationship between pulmonary microbes, immunity, and the lung epithelium is needed to improve outcomes. In this multicenter study, we collected 278 bronchoalveolar lavage (BAL) samples from 229 patients treated at 32 children's hospitals between 2014-2022. Using paired metatranscriptomes and human gene expression data, we identified 4 patient clusters with varying BAL composition. Among those requiring respiratory support prior to sampling, in-hospital mortality varied from 22-60% depending on the cluster (p=0.007). The most common patient subtype, Cluster 1, showed a moderate quantity and high diversity of commensal microbes with robust metabolic activity, low rates of infection, gene expression indicating alveolar macrophage predominance, and low mortality. The second most common cluster showed a very high burden of airway microbes, gene expression enriched for neutrophil signaling, frequent bacterial infections, and moderate mortality. Cluster 3 showed significant depletion of commensal microbes, a loss of biodiversity, gene expression indicative of fibroproliferative pathways, increased viral and fungal pathogens, and high mortality. Finally, Cluster 4 showed profound microbiome depletion with enrichment of Staphylococci and viruses, gene expression driven by lymphocyte activation and cellular injury, and the highest mortality. BAL clusters were modeled with a random forest classifier and reproduced in a geographically distinct validation cohort of 57 patients from The Netherlands, recapitulating similar cluster-based mortality differences (p=0.022). Degree of antibiotic exposure was strongly associated with depletion of BAL microbes and enrichment of fungi. Potential pathogens were parsed from all detected microbes by analyzing each BAL microbe relative to the overall microbiome composition, which yielded increased sensitivity for numerous previously occult pathogens. These findings support personalized interpretation of the pulmonary microenvironment in pediatric HCT, which may facilitate biology-targeted interventions to improve outcomes.

Evaluation of an established colorectal robotic programme at an NHS district general hospital: audit of outcomes and systematic review of published data.

INTRODUCTION: Robotic-assisted surgery (RAS) offers potential advantages over traditional surgical approaches. This study aimed to assess outcomes from a district general hospital (DGH) robotic colorectal programme against published data. MATERIALS AND METHODS: The robotic programme was established following simulator, dry/wet lab training, and proctoring. We performed a case series analysing technical, patient, and oncological outcomes extracted from a prospective database of colorectal RAS cases (2015-2022). A registered systematic review (PROSPERO CRD42022300773; PubMed, Web of Science, EMBASE) of single-centre colorectal series from established robotic centres (n>200 cases) was completed and compared to local data using descriptive summary statistics. Risk of bias assessment was performed using an adapted version of the Cochrane ROBINS-I tool. RESULTS: Two hundred thirty-two RAS cases were performed including 122 anterior resections, 56 APERs, 19 rectopexies, and 15 Hartmann's procedures. The median duration was 325 (IQR 265-400) min. Blood loss was < 100 ml in 97% of cases with 2 (0.9%) cases converted to open. Complications (Clavien-Dindo 3-5) occurred in 19 (8%) patients, with 3 (1.3%) deaths in < 30 days. Length of stay was 7 (IQR 5-11) days. In 169 rectal cancer cases, there were 9 (5.3%) cases with a positive circumferential or distal margin and lymph node yield of 17 (IQR 13-24). A systematic review of 1648 abstracts identified 13 studies from established robotic centres, totaling 4930 cases, with technical, patient, and oncological outcomes comparable to our own case series. CONCLUSIONS: Outcomes from our robotic colorectal programme at a UK DGH are comparable with the largest published case series from world-renowned centres. Training and proctoring together with rolling audit must accompany the expansion of robotic surgery to safeguard outcomes.

Population pharmacokinetic analyses for belzutifan to inform dosing considerations and labeling.

Belzutifan (Welireg, Merck & Co., Inc., Rahway, NJ, USA) is an oral, potent inhibitor of hypoxia-inducible factor 2α, approved for the treatment of certain patients with von Hippel-Lindau (VHL) disease-associated renal cell carcinoma (RCC), central nervous system hemangioblastomas, and pancreatic neuroendocrine tumors. It is primarily metabolized by the polymorphic uridine 5'-diphospho-glucuronosyltransferase (UGT) 2B17 and cytochrome (CYP) 2C19. A population pharmacokinetic (PK) model was built, using NONMEM version 7.3, based on demographics/PK data from three clinical pharmacology (food effect, formulation bridging, and genotype/race effect) and two clinical studies (phase I dose escalation/expansion in patients with RCC and other solid tumors; phase II in patients with VHL). Median (range) age for the combined studies was 55 years (19-84) and body weight was 73.6 kg (42.1-165.8). Belzutifan plasma PK was well-characterized by a linear two-compartment model with first-order absorption and elimination. For patients with VHL, the predicted geometric mean (% coefficient of variation) apparent clearance was 7.3 L/h (51%), apparent total volume of distribution was 130 L (35%), and half-life was 12.39 h (42%). There were no clinically relevant differences in belzutifan PK based on the individual covariates of age, sex, ethnicity, race, body weight, mild/moderate renal impairment, or mild hepatic impairment. In this model, dual UGT2B17 and CYP2C19 poor metabolizers (PMs) were estimated to have a 3.2-fold higher area under the plasma concentration-time curve compared to UGT2B17 extensive metabolizer and CYP2C19 non-PM patients. This population PK analysis enabled an integrated assessment of PK characteristics with covariate effects in the overall population and subpopulations for belzutifan labeling.

An Integrated Computational Analysis of High-Risk SNPs in Angiopoietin-like Proteins (ANGPTL3 and ANGPTL8) Reveals Perturbed Protein Dynamics Associated with Cancer.

Angiopoietin-like proteins (ANGPTL) constitute a family of eight proteins (1-8) which play a pivotal role in the regulation of various pathophysiological processes. The current study sought to identify high-risk, "non-synonymous, single-nucleotide polymorphisms" (nsSNPs) in both ANGPTL3 and ANGPTL8 to evaluate the role that these nsSNPs play in various types of cancer. We retrieved a total of 301 nsSNPs from various databases; 79 of these candidates constitute high-risk nsSNPs. Moreover, we identified eleven high-risk nsSNPs that cause various types of cancer: seven candidates for ANGPTL3 (L57H, F295L, L309F, K329M, R332L, S348C, and G409R) and four candidates for ANGPTL8 (P23L, R85W, R138S, and E148D). Protein-protein interaction analysis revealed a strong association of ANGPTL proteins with several tumor-suppressor proteins such as ITGB3, ITGAV, and RASSF5. 'Gene-expression profiling interactive analysis' (GEPIA) showed that expression of ANGPTL3 is significantly downregulated in five cancers: sarcoma (SARC); cholangio carcinoma (CHOL); kidney chromophobe carcinoma (KICH); kidney renal clear cell carcinoma (KIRC); and kidney renal papillary cell carcinoma (KIRP). GEPIA also showed that expression of ANGPTL8 remains downregulated in three cancers: CHOL; glioblastoma (GBM); and breast invasive carcinoma (BRCA). Survival rate analysis indicated that both upregulation and downregulation of ANGPTL3 and ANGPTL8 leads to low survival rates in various types of cancer. Overall, the current study revealed that both ANGPTL3 and ANGPTL8 constitute potential prognostic biomarkers for cancer; moreover, nsSNPs in these proteins might lead to the progression of cancer. However, further in vivo investigation will be helpful to validate the role of these proteins in the biology of cancer.

minSNPs: an R package for the derivation of resolution-optimised SNP sets from microbial genomic data.

Here, we present the R package, minSNPs. This is a re-development of a previously described Java application named Minimum SNPs. MinSNPs assembles resolution-optimised sets of single nucleotide polymorphisms (SNPs) from sequence alignments such as genome-wide orthologous SNP matrices. MinSNPs can derive sets of SNPs optimised for discriminating any user-defined combination of sequences from all others. Alternatively, SNP sets may be optimised to determine all sequences from all other sequences, i.e., to maximise diversity. MinSNPs encompasses functions that facilitate rapid and flexible SNP mining, and clear and comprehensive presentation of the results. The minSNPs' running time scales in a linear fashion with input data volume and the numbers of SNPs and SNPs sets specified in the output. MinSNPs was tested using a previously reported orthologous SNP matrix of Staphylococcus aureus and an orthologous SNP matrix of 3,279 genomes with 164,335 SNPs assembled from four S. aureus short read genomic data sets. MinSNPs was shown to be effective for deriving discriminatory SNP sets for potential surveillance targets and in identifying SNP sets optimised to discriminate isolates from different clonal complexes. MinSNPs was also tested with a large Plasmodium vivax orthologous SNP matrix. A set of five SNPs was derived that reliably indicated the country of origin within three south-east Asian countries. In summary, we report the capacity to assemble comprehensive SNP matrices that effectively capture microbial genomic diversity, and to rapidly and flexibly mine these entities for optimised marker sets.