A secondary mechanism of action for triazole antifungals in Aspergillus fumigatus mediated by hmg1

Triazole antifungals function as ergosterol biosynthesis inhibitors and are frontline therapy for invasive fungal infections, such as invasive aspergillosis. The primary mechanism of action of triazoles is through the specific inhibition of a cytochrome P450 14-α-sterol demethylase enzyme, Cyp51A/B, resulting in depletion of cellular ergosterol. Here, we uncover a clinically relevant secondary mechanism of action for triazoles within the ergosterol biosynthesis pathway. We provide evidence that triazole-mediated inhibition of Cyp51A/B activity generates sterol intermediate perturbations that are likely decoded by the sterol sensing functions of HMG-CoA reductase and Insulin-Induced Gene orthologs as increased pathway activity. This, in turn, results in negative feedback regulation of HMG-CoA reductase, the rate-limiting step of sterol biosynthesis. We also provide evidence that HMG-CoA reductase sterol sensing domain mutations previously identified as generating resistance in clinical isolates of Aspergillus fumigatus partially disrupt this triazole-induced feedback. Therefore, our data point to a secondary mechanism of action for the triazoles: induction of HMG-CoA reductase negative feedback for downregulation of ergosterol biosynthesis pathway activity. Abrogation of this feedback through acquired mutations in the HMG-CoA reductase sterol sensing domain diminishes triazole antifungal activity against fungal pathogens and underpins HMG-CoA reductase-mediated resistance.

TNIK inhibition sensitizes TNIK-overexpressing lung squamous cell carcinoma to radiotherapy

Most patients with lung squamous cell carcinoma (LSCC) undergo chemotherapy, radiotherapy, and adjuvant immunotherapy for locally advanced disease. The efficacy of these treatments is still limited due to dose-limiting toxicity or locoregional recurrence. New combination approaches and targets such as actionable oncogenic drivers are needed to advance treatment options for LSCC patients. Moreover, other options for chemotherapy-ineligible patients are also limited. As such there is a critical need for the development of selective and potent chemoradiosensitizers for locally advanced LSCC. Here, we investigated inhibiting TRAF2 and NCK-interacting protein kinase (TNIK), which is amplified in 40% of LSCC patients, as a strategy to sensitize LSCC tumors to chemo- and radiotherapy. Employing a range of human LSCC cell lines and the TNIK inhibitor NCB-0846, we investigated the potential of TNIK as a chemo- and radiosensitizing target with in vitro and in vivo preclinical models. The combination of NCB-0846 with cisplatin or etoposide was at best additive. Interestingly, pre-treating LSCC cells with NCB-0846 prior to ionizing radiation (IR) potentiated the cytotoxicity of IR in a TNIK-specific fashion. Characterization of the radiosensitization mechanism suggested that TNIK inhibition may impair the DNA damage response and promote mitotic catastrophe in irradiated cells. In a subcutaneous xenograft in vivo model, pretreatment with NCB-0846 significantly enhanced the efficacy of IR and caused elevated necrosis in TNIKhigh LK2 tumors but not TNIKlow KNS62 tumors. Overall, these results indicate that TNIK inhibition may be a promising strategy to increase the efficacy of radiotherapy in LSCC patients with high TNIK expression.

Prognostic Significance of Immune Cell Infiltration in Muscle-invasive Bladder Cancer Treated with Definitive Chemoradiation: A Secondary Analysis of RTOG 0524 and RTOG 0712

Chemoradiation therapy (CRT) is a treatment for muscle-invasive bladder cancer (MIBC). Using a novel transcriptomic profiling panel, we validated prognostic immune biomarkers to CRT using 70 pretreatment tumor samples from prospective trials of MIBC (NRG/RTOG 0524 and 0712). Disease-free survival (DFS) and overall survival (OS) were estimated via the Kaplan-Meier method and stratified by genes correlated with immune cell activation. Cox proportional-hazards models were used to assess group differences. Clustering of gene expression profiles revealed that the cluster with high immune cell content was associated with longer DFS (hazard ratio [HR] 0.53, 95% confidence interval [CI] 0.26-1.10; p = 0.071) and OS (HR 0.48, 95% CI 0.24-0.97; p = 0.040) than the cluster with low immune cell content. Higher expression of T-cell infiltration genes (CD8A and ICOS) was associated with longer DFS (HR 0.40, 95% CI 0.21-0.75; p = 0.005) and OS (HR 0.49, 95% CI 0.25-0.94; p = 0.033). Higher IDO1 expression (IFNγ signature) was also associated with longer DFS (HR 0.44, 95% CI 0.24-0.88; p = 0.021) and OS (HR 0.49, 95% CI 0.24-0.99; p = 0.048). These findings should be validated in prospective CRT trials that include biomarkers, particularly for trials incorporating immunotherapy for MIBC. PATIENT SUMMARY: We analyzed patient samples from two clinical trials (NRG/RTOG 0524 and 0712) of chemoradiation for muscle-invasive bladder cancer using a novel method to assess immune cells in the tumor microenvironment. Higher expression of genes associated with immune activation and high overall immune-cell content were associated with better disease-free survival and overall survival for patients treated with chemoradiation.

Dual targeting of glutamine and serine metabolism in acute myeloid leukemia

Acute myeloid leukemia (AML) is a heterogeneous hematological malignancy characterized by disrupted blood cell production and function. Recent investigations have highlighted the potential of targeting glutamine metabolism as a promising therapeutic approach for AML. Asparaginases, enzymes that deplete circulating glutamine and asparagine, are approved for the treatment of acute lymphoblastic leukemia, but are also under investigation in AML, with promising results. We previously reported an elevation in plasma serine levels following treatment with Erwinia-derived asparaginase (also called crisantaspase). This led us to hypothesize that AML cells initiate the de novo serine biosynthesis pathway in response to crisantaspase treatment and that inhibiting this pathway in combination with crisantaspase would enhance AML cell death. Here we report that in AML cell lines, treatment with the clinically available crisantaspase, Rylaze, upregulates the serine biosynthesis enzymes phosphoglycerate dehydrogenase (PHGDH) and phosphoserine aminotransferase (PSAT1) through activation of the Amino Acid Response (AAR) pathway, a cellular stress response mechanism that regulates amino acid metabolism and protein synthesis under conditions of nutrient limitation. Inhibition of serine biosynthesis through CRISPR-Cas9-mediated knockout of PHGDH resulted in a ~250-fold reduction in the half-maximal inhibitory concentration (IC50) for Rylaze, indicating heightened sensitivity to crisantaspase therapy. Treatment of AML cells with a combination of Rylaze and a small molecule inhibitor of PHGDH (BI4916) revealed synergistic anti-proliferative effects in both cell lines and primary AML patient samples. Rylaze-BI4916 treatment in AML cell lines led to the inhibition of cap-dependent mRNA translation and protein synthesis, as well as a marked decrease in intracellular glutathione levels, a critical cellular antioxidant. Collectively, our results highlight the clinical potential of targeting serine biosynthesis in combination with crisantaspase as a novel therapeutic strategy for AML.

Identification of a Predictive Genomic Biomarker for Prostate-directed Therapy in Synchronous Low-volume Metastatic Castration-sensitive Prostate Cancer

BACKGROUND

Standard of care management for synchronous metastatic castration-sensitive prostate cancer (mCSPC) includes androgen deprivation therapy with a second-generation antiandrogen therapy and/or docetaxel. Recently, randomized data have demonstrated that prostate-directed therapy (PDT) is associated with an improvement in overall survival (OS) among patients with low-volume metastatic disease. Tumor genomics represents an additional dimension to define the clinical trajectory of patients with mCSPC.

OBJECTIVE

To evaluate a high-risk (HiRi) genomic signature to predict the benefit from PDT.

DESIGN, SETTING, AND PARTICIPANTS

We performed a single-institution retrospective review of men with synchronous low-volume mCSPC who underwent DNA panel sequencing of their tumor. Patients were classified according to the presence of HiRi mutation including pathogenic mutations in TP53, ATM, BRCA1, BRCA2, or Rb1.

OUTCOME MEASUREMENTS AND STATISTICAL ANALYSIS

The primary endpoint was to determine the effect of PDT on OS in patients with and without a HiRi mutation. A survival analysis was performed with the Kaplan-Meier method compared with log-rank test and multivariable Cox regression. The interaction between HiRi mutation and PDT was evaluated.

RESULTS AND LIMITATIONS

A total of 101 patients with synchronous low-volume CSPC were included with a median follow-up of 44 mo. Approximately half of patients were found to have a HiRi pathogenic mutation (49%). Patients with HiRi mutations demonstrated median OS of 73 versus 66.8 mo (p = 0.3) for no PDT versus PDT. Conversely, patients without a HiRi mutation demonstrated a significant improvement in OS of 60 versus 105.3 mo (p < 0.001) for no PDT versus PDT. The p value for interaction for OS between PDT and HiRi mutation was statistically significant (p < 0.001). Limitations include the retrospective nature of the study.

CONCLUSIONS

Here, we have identified a HiRi genomic biomarker that appears predictive for the lack of benefit from PDT in men with synchronous low-volume mCSPC. Further work validating these results is warranted.

PATIENT SUMMARY

In this report, we evaluated a high-risk genomic biomarker to predict the benefit from prostate-directed therapy for men with synchronous low-volume metastatic castration-sensitive prostate cancer. We found that men without a high-risk mutation appear to experience a greater clinical benefit from prostate-directed therapy than those with a high-risk mutation.

TP53 structure-function relationships in metastatic castrate-sensitive prostate cancer and the impact of APR-246 treatment

PURPOSE

Despite well-informed work in several malignancies, the phenotypic effects of TP53 mutations in metastatic castration-sensitive prostate cancer (mCSPC) progression and metastasis are not clear. We characterized the structure-function and clinical impact of TP53 mutations in mCSPC.

PATIENTS AND METHODS

We performed an international retrospective review of men with mCSPC who underwent next-generation sequencing and were stratified according to TP53 mutational status and metastatic burden. Clinical outcomes included radiographic progression-free survival (rPFS) and overall survival (OS) evaluated with Kaplan-Meier and multivariable Cox regression. We also utilized isogenic cancer cell lines to assess the effect of TP53 mutations and APR-246 treatment on migration, invasion, colony formation in vitro, and tumor growth in vivo. Preclinical experimental observations were compared using t-tests and ANOVA.

RESULTS

Dominant-negative (DN) TP53 mutations were enriched in patients with synchronous (vs. metachronous) (20.7% vs. 6.3%, p < 0.01) and polymetastatic (vs. oligometastatic) (14.4% vs. 7.9%, p < 0.01) disease. On multivariable analysis, DN mutations were associated with worse rPFS (hazards ratio [HR] = 1.97, 95% confidence interval [CI]: 1.31-2.98) and overall survival [OS] (HR = 2.05, 95% CI: 1.14-3.68) compared to TP53 wild type (WT). In vitro, 22Rv1 TP53 R175H cells possessed stronger migration, invasion, colony formation ability, and cellular movement pathway enrichment in RNA sequencing analysis compared to 22Rv1 TP53 WT cells. Treatment with APR-246 reversed the effects of TP53 mutations in vitro and inhibited 22Rv1 TP53 R175H tumor growth in vivo in a dosage-dependent manner.

CONCLUSIONS

DN TP53 mutations correlated with worse prognosis in prostate cancer patients and higher metastatic potential, which could be counteracted by APR-246 treatment suggesting a potential future therapeutic avenue.

Clinical and Genomic Differences Between Advanced Molecular Imaging-detected and Conventional Imaging-detected Metachronous Oligometastatic Castration-sensitive Prostate Cancer

In metastatic castration-sensitive prostate cancer (mCSPC), disease volume plays an integral role in guiding treatment recommendations, including selection of docetaxel therapy, metastasis-directed therapy, and radiation to the prostate. Although there are multiple definitions of disease volume, they have commonly been studied in the context of metastases detected via conventional imaging (CIM). One such numeric definition of disease volume, termed oligometastasis, is heavily dependent on the sensitivity of the imaging modality. We performed an international multi-institutional retrospective review of men with metachronous oligometastatic CSPC (omCSPC), detected via either advanced molecular imaging alone (AMIM) or CIM. Patients were compared with respect to clinical and genomic features using the Mann-Whitney U test, Pearson's χ2 test, and Kaplan-Meier overall survival (OS) analyses with a log-rank test. A total of 295 patients were included for analysis. Patients with CIM-omCSPC had significantly higher Gleason grade group (p = 0.032), higher prostate-specific antigen at omCSPC diagnosis (8.0 vs 1.7 ng/ml; p < 0.001), more frequent pathogenic TP53 mutations (28% vs 17%; p = 0.030), and worse 10-yr OS (85% vs 100%; p < 0.001). This is the first report of clinical and biological differences between AMIM-detected and CIM-detected omCSPC. Our findings are particularly important for ongoing and planned clinical trials in omCSPC. PATIENT SUMMARY: Metastatic prostate cancer with just a few metastases only detected via newer scanning methods (called molecular imaging) is associated with fewer high-risk DNA mutations and better survival in comparison to metastatic cancer detected via conventional scan methods.

Cochlear organoids reveal transcriptional programs of postnatal hair cell differentiation from supporting cells

We explore the changes in chromatin accessibility and transcriptional programs for cochlear hair cell differentiation from postmitotic supporting cells using organoids from postnatal cochlea. The organoids contain cells with transcriptional signatures of differentiating vestibular and cochlear hair cells. Construction of trajectories identifies Lgr5+ cells as progenitors for hair cells, and the genomic data reveal gene regulatory networks leading to hair cells. We validate these networks, demonstrating dynamic changes both in expression and predicted binding sites of transcription factors (TFs) during organoid differentiation. We identify known regulators of hair cell development, Atoh1, Pou4f3, and Gfi1, and the analysis predicts the regulatory factors Tcf4, an E-protein and heterodimerization partner of Atoh1, and Ddit3, a CCAAT/enhancer-binding protein (C/EBP) that represses Hes1 and activates transcription of Wnt-signaling-related genes. Deciphering the signals for hair cell regeneration from mammalian cochlear supporting cells reveals candidates for hair cell (HC) regeneration, which is limited in the adult.

Single-cell RNA sequencing reveals peripheral blood mononuclear immune cell landscape associated with operational tolerance in a kidney transplant recipient

Operational tolerance (OT) after kidney transplantation is defined as stable graft acceptance without the need for immunosuppression therapy. However, it is not clear which cellular and molecular pathways are driving tolerance in these patients. In this first-of-its-kind pilot study, we assessed the immune landscape associated with OT using single-cell analyses. Peripheral mononuclear cells from a kidney transplant recipient with OT (Tol), 2 healthy individuals (HC), and a kidney transplant recipient with normal kidney function on standard-of-care immunosuppression (SOC) were evaluated. The immune landscape of the Tol was drastically different from that of SOC and emerged closer to the profile of HC. TCL1A+ naive B cells and LSGAL1+ regulatory T cells (Tregs) were in higher proportions in Tol. We were unable to identify the Treg subcluster in SOC. The ligand-receptor analysis in HC and Tol identified interactions between B cells, and Tregs that enhance the proliferation and suppressive function of Tregs. SOC reported the highest proportion of activated B cells with more cells in the G2M phase. Our single-cell RNA sequencing study identified the mediators of tolerance; however, it emphasizes the requirement of similar investigations on a larger cohort to reaffirm the role of immune cells in tolerance.

Novel Host Pathways Govern Epithelial Cell Invasion of Aspergillus fumigatus

Invasive aspergillosis is initiated when Aspergillus fumigatus adheres to and invades the pulmonary epithelial cells that line the airways and alveoli. To gain deeper insight into how pulmonary epithelial cells respond to A. fumigatus invasion, we used transcriptome sequencing (RNA-seq) to determine the transcriptional response of the A549 type II alveolar epithelial cell line to infection with strains CEA10 and Af293, two clinical isolates of A. fumigatus. Upstream regulator analysis of the data indicated that while both strains activated virtually identical host cell signaling pathways after 16 h of infection, only strain CEA10 activated these pathways after 6 h of infection. Many of the pathways that were predicted to be activated by A. fumigatus, including the tumor necrosis factor (TNF), interleukin-1α (IL-1α), IL-1β, IL-17A, Toll-like receptor 2 (TLR2), and TLR4 pathways, are known to be critical for the host defense against this fungus. We also found that the platelet-derived growth factor BB (PDGF BB) and progesterone receptor (PGR) pathways were activated by A. fumigatus. Using pharmacologic inhibitors, we determined that blocking the PDGF receptor or PGR inhibited the endocytosis of both strains of A. fumigatus in an additive manner. Both the PDGF BB and PGR pathways are also predicted to be activated by infection of A549 cells with other molds, such as Rhizopus delemar and Rhizopus oryzae. Thus, these pathways may represent a common response of pulmonary epithelial cells to mold infection. IMPORTANCE Invasive aspergillosis is a deadly invasive fungal infection that initiates when Aspergillus fumigatus spores are inhaled and come into contact with the epithelial cells that line the airways and alveoli. Understanding this fungus-host interaction is important for the development of novel therapeutics. To gain a deeper understanding of how these airway epithelial cells respond to A. fumigatus during infection, we used RNA-seq to determine the transcriptional response of alveolar epithelial cells to infection with two different clinical isolates of A. fumigatus. Our analysis identified new host response pathways that have not previously been tied to infection with A. fumigatus. Pharmacological inhibition of two of these pathways inhibited the ability of A. fumigatus to invade airway epithelial cells. These two pathways are also predicted to be activated by infection with other filamentous fungi. Thus, these pathways may represent a common response of alveolar epithelial cells to mold infection.

Proof of concept: Detection of cell free RNA from EDTA plasma in patients with lung cancer and non-cancer patients

INTRODUCTION

Nucleic acid sequencing technologies have advanced significantly in recent years, thereby allowing for the development of liquid biopsies as new means to detect cancer biomarkers and cancer heterogenicity. Most of the assays available, clinically, focus on cell free DNA (cfDNA), however, cell free RNA (cfRNA) is also present. cfRNA has the potential to complement and improve cancer detection especially in cancers like lung cancer, which are usually only diagnosed at late stages and therefore have poor long-term survival outcomes.

METHODS

Remnant EDTA plasma was collected from lung cancer patients and non-cancer individuals at the University of Maryland Medical Center. RNA was extracted and processed for next generation sequencing with a tagmentation-based library preparation approach.

RESULTS

cfRNA was successfully extracted and sequenced from 52 EDTA-treated plasma samples with volumes as low as 1.5 mL. This quantity was sufficient to prepare libraries with the length of libraries averaging from 264 bp to 381 bp and resulted in over 2.2 to 3.6 million total sequence reads respectively. Sequential dilution of cfRNA samples from healthy individuals indicated that the starting cfRNA concentration influenced the detection of differentially expressed genes.

CONCLUSIONS

This proof-of-concept study provides a framework for screening cfRNA for identifying biomarkers for early detection of lung cancer (and other cancers), using minimal amounts of samples (1.5 mL) from standard EDTA 3-mL collection tubes routinely used for patient care. Further studies in large populations are required to establish limit of detection and other parameters including precision, accuracy, sensitivity, and specificity, to standardize this method.

Transcriptomic and clinical heterogeneity of metastatic disease timing within metastatic castration-sensitive prostate cancer

BACKGROUND

Metastatic castration-sensitive prostate cancer (mCSPC) is commonly classified into high- and low-volume subgroups which have demonstrated differential biology, prognosis, and response to therapy. Timing of metastasis has similarly demonstrated differences in clinical outcomes; however, less is known about any underlying biologic differences between these disease states. Herein, we aim to compare transcriptomic differences between synchronous and metachronous mCSPC and identify any differential responses to therapy.

PATIENTS AND METHODS

We performed an international multi-institutional retrospective review of men with mCSPC who completed RNA expression profiling evaluation of their primary tumor. Patients were stratified according to disease timing (synchronous versus metachronous). The primary endpoint was to identify differences in transcriptomic profiles between disease timing. The median transcriptomic scores between groups were compared with the Mann-Whitney U test. Secondary analyses included determining clinical and transcriptomic variables associated with overall survival (OS) from the time of metastasis. Survival analysis was carried out with the Kaplan-Meier method and multivariable Cox regression.

RESULTS

A total of 252 patients were included with a median follow-up of 39.6 months. Patients with synchronous disease experienced worse 5-year OS (39% versus 79%; P < 0.01) and demonstrated lower median androgen receptor (AR) activity (11.78 versus 12.64; P < 0.01) and hallmark androgen response (HAR; 3.15 versus 3.32; P < 0.01). Multivariable Cox regression identified only high-volume disease [hazard ratio (HR) = 4.97, 95% confidence interval (CI) 2.71-9.10; P < 0.01] and HAR score (HR = 0.51, 95% CI 0.28-0.88; P = 0.02) significantly associated with OS. Finally, patients with synchronous (HR = 0.47, 95% CI 0.30-0.72; P < 0.01) but not metachronous (HR = 1.37, 95% CI 0.50-3.92; P = 0.56) disease were found to have better OS with AR and non-AR combination therapy as compared with monotherapy (P value for interaction = 0.05).

CONCLUSIONS

We have demonstrated a potential biologic difference between metastatic timing of mCSPC. Specifically, for patients with low-volume disease, those with metachronous low-volume disease have a more hormone-dependent transcriptional profile and exhibit a better prognosis than synchronous low-volume disease.

WNT Pathway Mutations in Metachronous Oligometastatic Castration-Sensitive Prostate Cancer

PURPOSE

WNT signaling is a cellular pathway that has been implicated in the development and progression of prostate cancer. Oligometastatic castration-sensitive prostate cancer (omCSPC) represents a unique state of disease in which metastasis-directed therapy (MDT) has demonstrated improvement in progression-free survival. Herein, we investigate the clinical implications of genomic alterations in the WNT signaling cascade in men with omCSPC.

METHODS AND MATERIALS

We performed an international multi-institutional retrospective study of 277 men with metachronous omCSPC who underwent targeted DNA sequencing of their primary/metastatic tumor. Patients were classified by presence or absence of pathogenic WNT pathway mutations (in the genes APC, RNF43, and CTNNB1). Pearson χ2 and Mann-Whitney U tests were used to determine differences in clinical factors between genomic strata. Kaplan-Meier survival curves were generated for radiographic progression-free survival and overall survival, stratified according to WNT pathway mutation status.

RESULTS

A pathogenic WNT pathway mutation was detected in 11.2% of patients. Patients with WNT pathway mutations were more likely to have visceral metastases (22.6% vs 2.8%; P < .01) and less likely to have regional lymph node metastases (29.0% vs 50.4%; P = .02). At time of oligometastasis, these patients were treated with MDT alone (33.9%), MDT + limited course of systemic therapy (20.6%), systemic therapy alone (22.4%), or observation (defined as no treatment for ≥6 months after metastatic diagnosis). Multivariable cox regression demonstrated WNT pathway mutations associated with significantly worse overall survival (hazard ratio, 3.87; 95% confidence interval, 1.25-12.00).

CONCLUSIONS

Somatic WNT pathway alterations are present in approximately 11% of patients with omCSPC and are associated with an increased likelihood of visceral metastases. Although these patients have a worse natural history, they may benefit from MDT.

Comparative transcriptomics from intestinal cells of permissive and non-permissive hosts during Ancylostoma ceylanicum infection reveals unique signatures of protection and host specificity

Soil-transmitted nematodes (STNs) place a tremendous burden on health and economics worldwide with an estimate of at least 1.5 billion people, or 24% of the population, being infected with at least 1 STN globally. Children and pregnant women carry the heavier pathological burden, and disease caused by the blood-feeding worm in the intestine can result in anaemia and delays in physical and intellectual development. These parasites are capable of infecting and reproducing in various host species, but what determines host specificity remains unanswered. Identifying the molecular determinants of host specificity would provide a crucial breakthrough towards understanding the biology of parasitism and could provide attractive targets for intervention. To investigate specificity mechanisms, members of the hookworm genus Ancylostoma provide a powerful system as they range from strict specialists to generalists. Using transcriptomics, differentially expressed genes (DEGs) in permissive (hamster) and non-permissive (mouse) hosts at different early time points during infection with A. ceylanicum were examined. Analysis of the data has identified unique immune responses in mice, as well as potential permissive signals in hamsters. Specifically, immune pathways associated with resistance to infection are upregulated in the non-permissive host, providing a possible protection mechanism that is absent in the permissive host. Furthermore, unique signatures of host specificity that may inform the parasite that it has invaded a permissive host were identified. These data provide novel insight into the tissue-specific gene expression differences between permissive and non-permissive hosts in response to hookworm infection.

Single nuclei transcriptomics delineates complex immune and kidney cell interactions contributing to kidney allograft fibrosis

Chronic allograft dysfunction (CAD), characterized histologically by interstitial fibrosis and tubular atrophy, is the major cause of kidney allograft loss. Here, using single nuclei RNA sequencing and transcriptome analysis, we identified the origin, functional heterogeneity, and regulation of fibrosis-forming cells in kidney allografts with CAD. A robust technique was used to isolate individual nuclei from kidney allograft biopsies and successfully profiled 23,980 nuclei from five kidney transplant recipients with CAD and 17,913 nuclei from three patients with normal allograft function. Our analysis revealed two distinct states of fibrosis in CAD; low and high extracellular matrix (ECM) with distinct kidney cell subclusters, immune cell types, and transcriptional profiles. Imaging mass cytometry analysis confirmed increased ECM deposition at the protein level. Proximal tubular cells transitioned to an injured mixed tubular (MT1) phenotype comprised of activated fibroblasts and myofibroblast markers, generated provisional ECM which recruited inflammatory cells, and served as the main driver of fibrosis. MT1 cells in the high ECM state achieved replicative repair evidenced by dedifferentiation and nephrogenic transcriptional signatures. MT1 in the low ECM state showed decreased apoptosis, decreased cycling tubular cells, and severe metabolic dysfunction, limiting the potential for repair. Activated B, T and plasma cells were increased in the high ECM state, while macrophage subtypes were increased in the low ECM state. Intercellular communication between kidney parenchymal cells and donor-derived macrophages, detected several years post-transplantation, played a key role in injury propagation. Thus, our study identified novel molecular targets for interventions aimed to ameliorate or prevent allograft fibrogenesis in kidney transplant recipients.

An optimized protocol for single nuclei isolation from clinical biopsies for RNA-seq

Single nuclei RNA sequencing (snRNA-seq) has evolved as a powerful tool to study complex human diseases. Single cell resolution enables the study of novel cell types, biological processes, cell trajectories, and cell–cell signaling pathways. snRNA-seq largely relies on the dissociation of intact nuclei from human tissues. However, the study of complex tissues using small core biopsies presents many technical challenges. Here, an optimized protocol for single nuclei isolation is presented for frozen and RNAlater preserved human kidney biopsies. The described protocol is fast, low cost, and time effective due to the elimination of cell sorting and ultra-centrifugation. Samples can be processed in 90 min or less. This method is effective for obtaining normal nuclei morphology without signs of structural damage. Using snRNA-seq, 16 distinct kidney cell clusters were recovered from normal and peri-transplant acute kidney injury allograft samples, including immune cell clusters. Quality control measurements demonstrated that these optimizations eliminated cellular debris and allowed for a high yield of high-quality nuclei and RNA for library preparation and sequencing. Cellular disassociation did not induce cellular stress responses, which recapitulated transcriptional patterns associated with standardized methods of nuclei isolation. Future applications of this protocol will allow for thorough investigations of small biobank biopsies, identifying cell-specific injury pathways and driving the discovery of novel diagnostics and therapeutic targets.

Translatome changes in acute myeloid leukemia cells post-exposure to pegcrisantaspase and venetoclax

The clinical outcomes of patients with AML treated with available therapy remains unsatisfactory. We recently reported that the BCL-2 inhibitor, venetoclax, synergized with pegcrisantaspase (PegC) and demonstrated remarkable in vivo efficacy in a preclinical model of AML with complex karyotype. Ven-PegC combination blocks synthesis of proteins in AML cells by inhibiting cap-dependent translation of mRNA. To further explore the impact of Ven-PegC on protein translation, we used polysome profiling and high-throughput RNA-seq to characterize Ven-PegC dependent changes to the translatome. Here we report that the translation of five mRNAs, including two microRNAs, one rRNA, and two mitochondrial genes was altered after exposure to all three treatments (Ven, PegC and Ven-PegC). We focused our translatome validation studies on six additional genes related to translational efficiency that were modified by Ven-PegC. Notably, Ven-PegC treatment increased the RNA translation and protein level of Tribbles homolog 3 (TRIB3), eukaryotic translation initiation factor 3 subunit C (eIF3C), doublesex and mab-3 related transcription factor 1 (DMRT1), salt inducible kinase 1 (SIK1). We validated the observed changes in gene/protein expression in vitro and confirmed our cell line-based studies in the bone marrow of an AML PDX model after Ven-PegC treatment. These results support examining alterations in the translatome post-chemotherapy to offer insight into drug mechanism of action and to inform future therapeutic decisions.

Tobacco Hornworm (Manduca sexta) caterpillars as a novel host model for the study of fungal virulence and drug efficacy

The two leading yeast pathogens of humans, Candida albicans and Cryptococcus neoformans, cause systemic infections in >1.4 million patients worldwide with mortality rates approaching 75%. It is thus imperative to study fungal virulence mechanisms, efficacy of antifungal drugs, and host response pathways. While this is commonly done in mammalian models, which are afflicted by ethical and practical concerns, invertebrate models, such as wax moth larvae and nematodes have been introduced over the last two decades. To complement existing invertebrate host models, we developed fifth instar caterpillars of the Tobacco Hornworm moth Manduca sexta as a novel host model. These caterpillars can be maintained at 37°C, are suitable for injections with defined amounts of yeast cells, and are susceptible to the most threatening yeast pathogens, including C. albicans, C. neoformans, C. auris, and C. glabrata. Importantly, fungal burden can be assessed daily throughout the course of infection in a single caterpillar’s feces and hemolymph. Infected caterpillars can be rescued by treatment with antifungal drugs. Notably, these animals are large enough for weight to provide a reliable and reproducible measure of fungal disease and to facilitate host tissue-specific expression analyses. M. sexta caterpillars combine a suite of parameters that make them suitable for the study of fungal virulence.

Venetoclax and pegcrisantaspase for complex karyotype acute myeloid leukemia

Complex karyotype acute myeloid leukemia (CK-AML) has a dismal outcome with current treatments, underscoring the need for new therapies. Here, we report synergistic anti-leukemic activity of the BCL-2 inhibitor venetoclax (Ven) and the asparaginase formulation Pegylated Crisantaspase (PegC) in CK-AML in vitro and in vivo. Ven-PegC combination inhibited growth of multiple AML cell lines and patient-derived primary CK-AML cells in vitro. In vivo, Ven-PegC showed potent reduction of leukemia burden and improved survival, compared with each agent alone, in a primary patient-derived CK-AML xenograft. Superiority of Ven-PegC, compared to single drugs, and, importantly, the clinically utilized Ven-azacitidine combination, was also demonstrated in vivo in CK-AML. We hypothesized that PegC-mediated plasma glutamine depletion inhibits 4EBP1 phosphorylation, decreases the expression of proteins such as MCL-1, whose translation is cap dependent, synergizing with the BCL-2 inhibitor Ven. Ven-PegC treatment decreased cellular MCL-1 protein levels in vitro by enhancing eIF4E-4EBP1 interaction on the cap-binding complex via glutamine depletion. In vivo, Ven-PegC treatment completely depleted plasma glutamine and asparagine and inhibited mRNA translation and cellular protein synthesis. Since this novel mechanistically-rationalized regimen combines two drugs already in use in acute leukemia treatment, we plan a clinical trial of the Ven-PegC combination in relapsed/refractory CK-AML.

Tissue Damage in Radiation-Induced Oral Mucositis Is Mitigated by IL-17 Receptor Signaling

Oral mucositis (OM) is a treatment-limiting adverse side effect of radiation and chemotherapy. Approximately 80% of patients undergoing radiotherapy (RT) for head and neck cancers (HNC) develop OM, representing a major unmet medical condition. Our understanding of the immunopathogenesis of OM is limited, due in part to the surprising paucity of information regarding healing mechanisms in the oral mucosa. RNAseq of oral tissue in a murine model that closely mimics human OM, showed elevated expression of IL-17 and related immune pathways in response to head and neck irradiation (HNI). Strikingly, mice lacking the IL-17 receptor (IL-17RA) exhibited markedly more severe OM. Restoration of the oral mucosa was compromised in Il17ra-/- mice and components associated with healing, including matrix metalloproteinase 3, 10 and IL-24 were diminished. IL-17 is typically associated with recruitment of neutrophils to mucosal sites following oral infections. Unexpectedly, in OM the absence of IL-17RA resulted in excessive neutrophil recruitment and immunopathology. Instead, neutrophil activation was IL-1R-driven in Il17ra-/- mice. Blockade of IL-1R and depletion of neutrophils lessened the severity of damage in these mice. Overall, we show IL-17 is protective in OM through multiple mechanisms including restoration of the damaged epithelia and control of the neutrophil response. We also present a clinically relevant murine model of human OM to improve mechanistic understanding and develop rational translational therapeutics.

Oral mucosal IL-22/STAT3 signaling licenses IL-17-mediated immunity to oral candidiasis

Oropharyngeal candidiasis (OPC) is an opportunistic infection of the oral mucosa caused by the commensal fungus Candida albicans. IL-17 and IL-22 both mediate antifungal immunity yet activate distinct downstream signaling pathways. While much is known about IL-17-dependent immunity in OPC, the activities of IL-22 are less well delineated. We show that induction of Il22 is independent of Dectin-1, CARD9 and aryl hydrocarbon receptor (AhR) and is driven by IL-23 and the C. albicans pore forming peptide candidalysin. Despite similar induction requirements and cellular sources, IL-22 and IL-17 function non-redundantly during OPC and exert opposing roles in neutrophil recruitment. The IL-22 and IL-17 receptors are required in anatomically distinct locations; loss of IL-22RA1 in the oral basal epithelial layer (BEL) but not the suprabasal epithelial layer (SEL) causes susceptibility to OPC, whereas IL-17RA is needed in the SEL. Our data reveal that IL-22 is a major activator of STAT3 in the BEL during OPC. Moreover, loss of STAT3 in the BEL but not the SEL renders mice susceptible to OPC. Transcriptional profiling of RNASeq data linked IL-22/STAT3 to oral epithelial cell proliferation and survival, but also, unexpectedly, to driving an IL-17 gene signature. We show that IL-22 acts on the BEL to replenish the IL-17RA-expressing SEL, thereby restoring the ability of the oral epithelium to respond to IL-17. Consequently, IL-22 signaling in BEL ‘licenses’ IL-17R signaling in the oral epithelium, revealing spatially distinct yet cooperative activities of IL-22 and IL-17 in oral candidiasis. This work also suggests that oral thrush in Jobs’ syndrome patients may be caused by STAT3 impairments in the oral epithelium, not just Th17 cells.

Best practices on the differential expression analysis of multi-species RNA-seq

Advances in transcriptome sequencing allow for simultaneous interrogation of differentially expressed genes from multiple species originating from a single RNA sample, termed dual or multi-species transcriptomics. Compared to single-species differential expression analysis, the design of multi-species differential expression experiments must account for the relative abundances of each organism of interest within the sample, often requiring enrichment methods and yielding differences in total read counts across samples. The analysis of multi-species transcriptomics datasets requires modifications to the alignment, quantification, and downstream analysis steps compared to the single-species analysis pipelines. We describe best practices for multi-species transcriptomics and differential gene expression.

Determining Aspergillus fumigatus transcription factor expression and function during invasion of the mammalian lung

To gain a better understanding of the transcriptional response of Aspergillus fumigatus during invasive pulmonary infection, we used a NanoString nCounter to assess the transcript levels of 467 A. fumigatus genes during growth in the lungs of immunosuppressed mice. These genes included ones known to respond to diverse environmental conditions and those encoding most transcription factors in the A. fumigatus genome. We found that invasive growth in vivo induces a unique transcriptional profile as the organism responds to nutrient limitation and attack by host phagocytes. This in vivo transcriptional response is largely mimicked by in vitro growth in Aspergillus minimal medium that is deficient in nitrogen, iron, and/or zinc. From the transcriptional profiling data, we selected 9 transcription factor genes that were either highly expressed or strongly up-regulated during in vivo growth. Deletion mutants were constructed for each of these genes and assessed for virulence in mice. Two transcription factor genes were found to be required for maximal virulence. One was rlmA, which is required for the organism to achieve maximal fungal burden in the lung. The other was sltA, which regulates of the expression of multiple secondary metabolite gene clusters and mycotoxin genes independently of laeA. Using deletion and overexpression mutants, we determined that the attenuated virulence of the ΔsltA mutant is due in part to decreased expression aspf1, which specifies a ribotoxin, but is not mediated by reduced expression of the fumigaclavine gene cluster or the fumagillin-pseruotin supercluster. Thus, in vivo transcriptional profiling focused on transcription factors genes provides a facile approach to identifying novel virulence regulators.

Although A. fumigatus causes the majority of cases of invasive aspergillosis, the function of most genes in its genome remains unknown. To identify genes encoding transcription factors that may be important for virulence, we used a NanoString nCounter to measure the mRNA levels of A. fumigatus transcription factor genes in the lungs of mice with invasive aspergillosis. The transcriptional profiling data indicate that the organism is exposed to nutrient limitation and stress during growth in the lungs, and that it responds by up-regulating genes that encode mycotoxins and secondary metabolites. In vitro, this response was most closely mimicked by growth in medium that was deficient in nitrogen, iron and/or zinc. Using the transcriptional profiling data, we identified two transcription factors that govern A. fumigatus virulence. These were RlmA, which is governs factors that enables the organism to proliferate maximally in the lung and SltA, which controls the production of mycotoxins and secondary metabolites.

Therapeutic implications of C. albicans-S. aureus mixed biofilm in a murine subcutaneous catheter model of polymicrobial infection

Biofilm-associated polymicrobial infections tend to be challenging to treat. Candida albicans and Staphylococcus aureus are leading pathogens due to their ability to form biofilms on medical devices. However, the therapeutic implications of their interactions in a host is largely unexplored. In this study, we used a mouse subcutaneous catheter model for in vivo-grown polymicrobial biofilms to validate our in vitro findings on C. albicans-mediated enhanced S. aureus tolerance to vancomycin in vivo. Comparative assessment of S. aureus recovery from catheters with single- or mixed-species infection demonstrated failure of vancomycin against S. aureus in mice with co-infected catheters. To provide some mechanistic insights, RNA-seq analysis was performed on catheter biofilms to delineate transcriptional modulations during polymicrobial infections. C. albicans induced the activation of the S. aureus biofilm formation network via down-regulation of the lrg operon, repressor of autolysis, and up-regulation of the ica operon and production of polysaccharide intercellular adhesin (PIA), indicating an increase in eDNA production, and extracellular polysaccharide matrix, respectively. Interestingly, virulence factors important for disseminated infections, and superantigen-like proteins were down-regulated during mixed-species infection, whereas capsular polysaccharide genes were up-regulated, signifying a strategy favoring survival, persistence and host immune evasion. In vitro follow-up experiments using DNA enzymatic digestion, lrg operon mutant strains, and confocal scanning microscopy confirmed the role of C. albicans-mediated enhanced eDNA production in mixed-biofilms on S. aureus tolerance to vancomycin. Combined, these findings provide mechanistic insights into the therapeutic implications of interspecies interactions, underscoring the need for novel strategies to overcome limitations of current therapies.

Sequencing of 53,831 diverse genomes from the NHLBI TOPMed Program

The Trans-Omics for Precision Medicine (TOPMed) programme seeks to elucidate the genetic architecture and biology of heart, lung, blood and sleep disorders, with the ultimate goal of improving diagnosis, treatment and prevention of these diseases. The initial phases of the programme focused on whole-genome sequencing of individuals with rich phenotypic data and diverse backgrounds. Here we describe the TOPMed goals and design as well as the available resources and early insights obtained from the sequence data. The resources include a variant browser, a genotype imputation server, and genomic and phenotypic data that are available through dbGaP (Database of Genotypes and Phenotypes)1. In the first 53,831 TOPMed samples, we detected more than 400 million single-nucleotide and insertion or deletion variants after alignment with the reference genome. Additional previously undescribed variants were detected through assembly of unmapped reads and customized analysis in highly variable loci. Among the more than 400 million detected variants, 97% have frequencies of less than 1% and 46% are singletons that are present in only one individual (53% among unrelated individuals). These rare variants provide insights into mutational processes and recent human evolutionary history. The extensive catalogue of genetic variation in TOPMed studies provides unique opportunities for exploring the contributions of rare and noncoding sequence variants to phenotypic variation. Furthermore, combining TOPMed haplotypes with modern imputation methods improves the power and reach of genome-wide association studies to include variants down to a frequency of approximately 0.01%.

Monocyte Subsets With High Osteoclastogenic Potential and Their Epigenetic Regulation Orchestrated by IRF8

Osteoclasts (OCs) are bone-resorbing cells formed by the serial fusion of monocytes. In mice and humans, three distinct subsets of monocytes exist; however, it is unclear if all of them exhibit osteoclastogenic potential. Here we show that in wild-type (WT) mice, Ly6C^hi and Ly6C^int monocytes are the primary source of OC formation when compared to Ly6C^- monocytes. Their osteoclastogenic potential is dictated by increased expression of signaling receptors and activation of preestablished transcripts, as well as de novo gain in enhancer activity and promoter changes. In the absence of interferon regulatory factor 8 (IRF8), a transcription factor important for myelopoiesis and osteoclastogenesis, all three monocyte subsets are programmed to display higher osteoclastogenic potential. Enhanced NFATc1 nuclear translocation and amplified transcriptomic and epigenetic changes initiated at early developmental stages direct the increased osteoclastogenesis in Irf8-deficient mice. Collectively, our study provides novel insights into the transcription factors and active cis-regulatory elements that regulate OC differentiation. © 2020 American Society for Bone and Mineral Research (ASBMR).

Cost effective, experimentally robust differential-expression analysis for human/mammalian, pathogen and dual-species transcriptomics

As sequencing read length has increased, researchers have quickly adopted longer reads for their experiments. Here, we examine 14 pathogen or host–pathogen differential gene expression data sets to assess whether using longer reads is warranted. A variety of data sets was used to assess what genomic attributes might affect the outcome of differential gene expression analysis including: gene density, operons, gene length, number of introns/exons and intron length. No genome attribute was found to influence the data in principal components analysis, hierarchical clustering with bootstrap support, or regression analyses of pairwise comparisons that were undertaken on the same reads, looking at all combinations of paired and unpaired reads trimmed to 36, 54, 72 and 101 bp. Read pairing had the greatest effect when there was little variation in the samples from different conditions or in their replicates (e.g. little differential gene expression). But overall, 54 and 72 bp reads were typically most similar. Given differences in costs and mapping percentages, we recommend 54 bp reads for organisms with no or few introns and 72 bp reads for all others. In a third of the data sets, read pairing had absolutely no effect, despite paired reads having twice as much data. Therefore, single-end reads seem robust for differential-expression analyses, but in eukaryotes paired-end reads are likely desired to analyse splice variants and should be preferred for data sets that are acquired with the intent to be community resources that might be used in secondary data analyses.

Repeated sampling facilitates within- and between-subject modeling of the human sperm transcriptome to identify dynamic and stress-responsive sncRNAs

Epidemiological studies from the last century have drawn strong associations between paternal life experiences and offspring health and disease outcomes. Recent studies have demonstrated sperm small non-coding RNA (sncRNA) populations vary in response to diverse paternal insults. However, for studies in retrospective or prospective human cohorts to identify changes in paternal germ cell epigenetics in association with offspring disease risk, a framework must first be built with insight into the expected biological variation inherent in human populations. In other words, how will we know what to look for if we don't first know what is stable and what is dynamic, and what is consistent within and between men over time? From sperm samples from a 'normative' cohort of healthy human subjects collected repeatedly from each subject over 6 months, 17 healthy male participants met inclusion criteria and completed donations and psychological evaluations of perceived stress monthly. sncRNAs (including miRNA, piRNA, and tRNA) isolated from mature sperm from these samples were subjected to Illumina small RNA sequencing, aligned to subtype-specific reference transcriptomes, and quantified. The repeated measures design allowed us to define both within- and between-subject variation in the expression of 254 miRNA, 194 tRNA, and 937 piRNA in sperm over time. We developed screening criteria to identify a subset of potential environmentally responsive 'dynamic' sperm sncRNA. Implementing complex modeling of the relationships between individual dynamic sncRNA and perceived stress states in these data, we identified 5 miRNA (including let-7f-5p and miR-181a-5p) and 4 tRNA that are responsive to the dynamics of prior stress experience and fit our established mouse model. In the current study, we aligned repeated sampling of human sperm sncRNA expression data with concurrent measures of perceived stress as a novel framework that can now be applied across a range of studies focused on diverse environmental factors able to influence germ cell programming and potentially impact offspring development.

Oral epithelial IL-22/STAT3 signaling licenses IL-17-mediated immunity to oral mucosal candidiasis

Oropharyngeal candidiasis (OPC; thrush) is an opportunistic infection caused by the commensal fungus Candida albicans Interleukin-17 (IL-17) and IL-22 are cytokines produced by type 17 lymphocytes. Both cytokines mediate antifungal immunity yet activate quite distinct downstream signaling pathways. While much is now understood about how IL-17 promotes immunity in OPC, the activities of IL-22 are far less well delineated. We show that, despite having similar requirements for induction from type 17 cells, IL-22 and IL-17 function nonredundantly during OPC. We find that the IL-22 and IL-17 receptors are required in anatomically distinct locations within the oral mucosa; loss of IL-22RA1 or signal transducer and activator of transcription 3 (STAT3) in the oral basal epithelial layer (BEL) causes susceptibility to OPC, whereas IL-17RA is needed in the suprabasal epithelial layer (SEL). Transcriptional profiling of the tongue linked IL-22/STAT3 not only to oral epithelial cell proliferation and survival but also, unexpectedly, to driving an IL-17-specific gene signature. We show that IL-22 mediates regenerative signals on the BEL that replenish the IL-17RA-expressing SEL, thereby restoring the ability of the oral epithelium to respond to IL-17 and thus to mediate antifungal events. Consequently, IL-22 signaling in BEL "licenses" IL-17 signaling in the oral mucosa, revealing spatially distinct yet cooperative activities of IL-22 and IL-17 in oral candidiasis.

Oral epithelial IL-22/STAT3 signaling licenses IL-17-mediated immunity to oral mucosal candidiasis

Oropharyngeal candidiasis (OPC) is an opportunistic infection of the oral mucosa caused by the commensal fungus Candida albicans. IL-17 and IL-22 both mediate antifungal immunity yet activate distinct downstream signaling pathways. While much is known about IL-17-dependent immunity in OPC, the activities of IL-22 are less well delineated. We show that induction of Il22 is independent of Dectin-1, CARD9 and aryl hydrocarbon receptor (AhR) and is driven by IL-23 and the C. albicans pore forming peptide candidalysin. Despite similar induction requirements and cellular sources, IL-22 and IL-17 function non-redundantly during OPC and exert opposing roles in neutrophil recruitment. The IL-22 and IL-17 receptors are required in anatomically distinct locations; loss of IL-22RA1 in the oral basal epithelial layer (BEL) but not the suprabasal epithelial layer (SEL) causes susceptibility to OPC, whereas IL-17RA is needed in the SEL. Our data reveal that IL-22 is a major activator of STAT3 in the BEL during OPC. Moreover, loss of STAT3 in the BEL but not the SEL renders mice susceptible to OPC. Transcriptional profiling of RNASeq data linked IL-22/STAT3 to oral epithelial cell proliferation and survival, but also, unexpectedly, to driving an IL-17 gene signature. We show that IL-22 acts on the BEL to replenish the IL-17RA-expressing SEL, thereby restoring the ability of the oral epithelium to respond to IL-17. Consequently, IL-22 signaling in BEL ‘licenses’ IL-17R signaling in the oral epithelium, revealing spatially distinct yet cooperative activities of IL-22 and IL-17 in oral candidiasis. This work also suggests that oral thrush in Jobs’ syndrome patients may be caused by STAT3 impairments in the oral epithelium, not just Th17 cells.

Evolutionary history of modern Samoans

Archaeological studies estimate the initial settlement of Samoa at 2,750 to 2,880 y ago and identify only limited settlement and human modification to the landscape until about 1,000 to 1,500 y ago. At this point, a complex history of migration is thought to have begun with the arrival of people sharing ancestry with Near Oceanic groups (i.e., Austronesian-speaking and Papuan-speaking groups), and was then followed by the arrival of non-Oceanic groups during European colonialism. However, the specifics of this peopling are not entirely clear from the archaeological and anthropological records, and is therefore a focus of continued debate. To shed additional light on the Samoan population history that this peopling reflects, we employ a population genetic approach to analyze 1,197 Samoan high-coverage whole genomes. We identify population splits between the major Samoan islands and detect asymmetrical gene flow to the capital city. We also find an extreme bottleneck until about 1,000 y ago, which is followed by distinct expansions across the islands and subsequent bottlenecks consistent with European colonization. These results provide for an increased understanding of Samoan population history and the dynamics that inform it, and also demonstrate how rapid demographic processes can shape modern genomes.

Detecting geospatial patterns of Plasmodium falciparum parasite migration in Cambodia using optimized estimated effective migration surfaces

BACKGROUND

Understanding the genetic structure of natural populations provides insight into the demographic and adaptive processes that have affected those populations. Such information, particularly when integrated with geospatial data, can have translational applications for a variety of fields, including public health. Estimated effective migration surfaces (EEMS) is an approach that allows visualization of the spatial patterns in genomic data to understand population structure and migration. In this study, we developed a workflow to optimize the resolution of spatial grids used to generate EEMS migration maps and applied this optimized workflow to estimate migration of Plasmodium falciparum in Cambodia and bordering regions of Thailand and Vietnam.

METHODS

The optimal density of EEMS grids was determined based on a new workflow created using density clustering to define genomic clusters and the spatial distance between genomic clusters. Topological skeletons were used to capture the spatial distribution for each genomic cluster and to determine the EEMS grid density; i.e., both genomic and spatial clustering were used to guide the optimization of EEMS grids. Model accuracy for migration estimates using the optimized workflow was tested and compared to grid resolutions selected without the optimized workflow. As a test case, the optimized workflow was applied to genomic data generated from P. falciparum sampled in Cambodia and bordering regions, and migration maps were compared to estimates of malaria endemicity, as well as geographic properties of the study area, as a means of validating observed migration patterns.

RESULTS

Optimized grids displayed both high model accuracy and reduced computing time compared to grid densities selected in an unguided manner. In addition, EEMS migration maps generated for P. falciparum using the optimized grid corresponded to estimates of malaria endemicity and geographic properties of the study region that might be expected to impact malaria parasite migration, supporting the validity of the observed migration patterns.

CONCLUSIONS

Optimized grids reduce spatial uncertainty in the EEMS contours that can result from user-defined parameters, such as the resolution of the spatial grid used in the model. This workflow will be useful to a broad range of EEMS users as it can be applied to analyses involving other organisms of interest and geographic areas.

Candida albicans quorum-sensing molecule farnesol modulates staphyloxanthin production and activates the thiol-based oxidative-stress response in Staphylococcus aureus

Microbial species utilize secreted-signaling molecules to coordinate their behavior. Our previous investigations demonstrated a key role for the Candida albicans-secreted quorum-sensing molecule farnesol in modulating Staphylococcus aureus response to antimicrobials in mixed biofilms. In this study, we aimed to provide mechanistic insights into the impact of farnesol on S. aureus within the context of inter-species interactions. To mimic biofilm dynamics, farnesol-sensitized S. aureus cells were generated via sequential farnesol exposure. The sensitized phenotype exhibited dramatic loss of the typical pigment, which we identified as staphyloxanthin, an important virulence factor synthesized by the Crt operon in S. aureus. Additionally, farnesol exposure exerted oxidative-stress as indicated by transcriptional analysis demonstrating alterations in redox-sensors and major virulence regulators. Paradoxically, the activated stress-response conferred S. aureus with enhanced tolerance to H₂O₂ and phagocytic killing. Since expression of enzymes in the staphyloxanthin biosynthesis pathway was not impacted by farnesol, we generated a theoretical-binding model which indicated that farnesol may block staphyloxanthin biosynthesis via competitive-binding to the CrtM enzyme crucial for staphyloxanthin synthesis, due to high structural similarity to the CrtM substrate. Finally, mixed growth with C. albicans was found to similarly induce S. aureus depigmentation, but not during growth with a farnesol-deficient C. albicans strain. Collectively, the findings demonstrate that a fungal molecule acts as a redox-cycler eliciting a bacterial stress response via activation of the thiol-based redox system under the control of global regulators. Therefore, farnesol-induced transcriptional modulations of key regulatory networks in S. aureus may modulate the pathogenesis of C. albicans-S. aureus co-infections.

Bioinformatic analysis of bacteria and host cell dual RNA-sequencing experiments

Bacterial pathogens subvert host cells by manipulating cellular pathways for survival and replication; in turn, host cells respond to the invading pathogen through cascading changes in gene expression. Deciphering these complex temporal and spatial dynamics to identify novel bacterial virulence factors or host response pathways is crucial for improved diagnostics and therapeutics. Dual RNA sequencing (dRNA-Seq) has recently been developed to simultaneously capture host and bacterial transcriptomes from an infected cell. This approach builds on the high sensitivity and resolution of RNA sequencing technology and is applicable to any bacteria that interact with eukaryotic cells, encompassing parasitic, commensal or mutualistic lifestyles. Several laboratory protocols have been presented that outline the collection, extraction and sequencing of total RNA for dRNA-Seq experiments, but there is relatively little guidance available for the detailed bioinformatic analyses required. This protocol outlines a typical dRNA-Seq experiment, based on a Chlamydia trachomatis-infected host cell, with a detailed description of the necessary bioinformatic analyses with currently available software tools.

Genomic structure and diversity of Plasmodium falciparum in Southeast Asia reveal recent parasite migration patterns

Estimates of Plasmodium falciparum migration may inform strategies for malaria elimination. Here we elucidate fine-scale parasite population structure and infer recent migration across Southeast Asia using identity-by-descent (IBD) approaches based on genome-wide single nucleotide polymorphisms called in 1722 samples from 54 districts. IBD estimates are consistent with isolation-by-distance. We observe greater sharing of larger IBD segments between artemisinin-resistant parasites versus sensitive parasites, which is consistent with the recent spread of drug resistance. Our IBD analyses reveal actionable patterns, including isolated parasite populations, which may be prioritized for malaria elimination, as well as asymmetrical migration identifying potential sources and sinks of migrating parasites.

Iron restriction inside macrophages regulates pulmonary host defense against Rhizopus species

Mucormycosis is a life-threatening respiratory fungal infection predominantly caused by Rhizopus species. Mucormycosis has incompletely understood pathogenesis, particularly how abnormalities in iron metabolism compromise immune responses. Here we show how, as opposed to other filamentous fungi, Rhizopus spp. establish intracellular persistence inside alveolar macrophages (AMs). Mechanistically, lack of intracellular swelling of Rhizopus conidia results in surface retention of melanin, which induces phagosome maturation arrest through inhibition of LC3-associated phagocytosis. Intracellular inhibition of Rhizopus is an important effector mechanism, as infection of immunocompetent mice with swollen conidia, which evade phagocytosis, results in acute lethality. Concordantly, AM depletion markedly increases susceptibility to mucormycosis. Host and pathogen transcriptomics, iron supplementation studies, and genetic manipulation of iron assimilation of fungal pathways demonstrate that iron restriction inside macrophages regulates immunity against Rhizopus. Our findings shed light on the pathogenetic mechanisms of mucormycosis and reveal the role of macrophage-mediated nutritional immunity against filamentous fungi.

Mucormycosis is a life-threatening respiratory fungal infection that typically occurs in patients with abnormalities in iron metabolism. Here the authors show that iron restriction inside the phagosome of macrophages is an essential component of the host defense against Rhizopus, the main species causing mucormycosis.

Alternative Splicing of Differentiated Myeloid Cell Transcripts after Infection by Anaplasma phagocytophilum Impacts a Selective Group of Cellular Programs

Eukaryotic proteome diversity exceeds that encoded within individual genes, and results in part from alternative splicing events of pre-messenger RNA. The diversity of these splicing events can shape the outcome in development and differentiation of normal tissues, and is important in pathogenic circumstances such as cancer and some heritable conditions. A role for alternative splicing of eukaryotic genes in response to viral and intracellular bacterial infections has only recently been recognized, and plays an important role in providing fitness for microbial survival, while potentially enhancing pathogenicity. Anaplasma phagocytophilum survives within mammalian neutrophils by reshaping transcriptional programs that govern cellular functions. We applied next generation RNAseq to ATRA-differentiated HL-60 cells established to possess transcriptional and functional responses similar to A. phagocytophilum-infected human neutrophils. This demonstrated an increase in transcripts with infection and high proportion of alternatively spliced transcript events (ASEs) for which predicted gene ontology processes were in part distinct from those identified by evaluation of single transcripts or gene-level analyses alone. The alternative isoforms are not on average shorter, and no alternative splicing in genes encoding spliceosome components is noted. Although not evident at gene-level analyses, individual spliceosome transcripts that impact nearly all spliceosome components were significantly upregulated. How the distinct GO processes predicted by ASEs are regulated by infection and whether they are relevant to fitness or pathogenicity of A. phagocytophilum should be addressed in more detailed studies.

Cross-species transcriptional analysis reveals conserved and host-specific neoplastic processes in mammalian glioma

Glioma is a unique neoplastic disease that develops exclusively in the central nervous system (CNS) and rarely metastasizes to other tissues. This feature strongly implicates the tumor-host CNS microenvironment in gliomagenesis and tumor progression. We investigated the differences and similarities in glioma biology as conveyed by transcriptomic patterns across four mammalian hosts: rats, mice, dogs, and humans. Given the inherent intra-tumoral molecular heterogeneity of human glioma, we focused this study on tumors with upregulation of the platelet-derived growth factor signaling axis, a common and early alteration in human gliomagenesis. The results reveal core neoplastic alterations in mammalian glioma, as well as unique contributions of the tumor host to neoplastic processes. Notable differences were observed in gene expression patterns as well as related biological pathways and cell populations known to mediate key elements of glioma biology, including angiogenesis, immune evasion, and brain invasion. These data provide new insights regarding mammalian models of human glioma, and how these insights and models relate to our current understanding of the human disease.

Association of a Novel Mutation in the Plasmodium falciparum Chloroquine Resistance Transporter With Decreased Piperaquine Sensitivity

Background

Amplified copy number in the plasmepsin II/III genes within Plasmodium falciparum has been associated with decreased sensitivity to piperaquine. To examine this association and test whether additional loci might also contribute, we performed a genome-wide association study of ex vivo P. falciparum susceptibility to piperaquine.

Methods

Plasmodium falciparum DNA from 183 samples collected primarily from Cambodia was genotyped at 33716 genome-wide single nucleotide polymorphisms (SNPs). Linear mixed models and random forests were used to estimate associations between parasite genotypes and piperaquine susceptibility. Candidate polymorphisms were evaluated for their association with dihydroartemisinin-piperaquine treatment outcomes in an independent dataset.

Results

Single nucleotide polymorphisms on multiple chromosomes were associated with piperaquine 90% inhibitory concentrations (IC90) in a genome-wide analysis. Fine-mapping of genomic regions implicated in genome-wide analyses identified multiple SNPs in linkage disequilibrium with each other that were significantly associated with piperaquine IC90, including a novel mutation within the gene encoding the P. falciparum chloroquine resistance transporter, PfCRT. This mutation (F145I) was associated with dihydroartemisinin-piperaquine treatment failure after adjusting for the presence of amplified plasmepsin II/III, which was also associated with decreased piperaquine sensitivity.

Conclusions

Our data suggest that, in addition to plasmepsin II/III copy number, other loci, including pfcrt, may also be involved in piperaquine resistance.

Streptococcus pneumoniae in the heart subvert the host response through biofilm-mediated resident macrophage killing

For over 130 years, invasive pneumococcal disease has been associated with the presence of extracellular planktonic pneumococci, i.e. diplococci or short chains in affected tissues. Herein, we show that Streptococcus pneumoniae that invade the myocardium instead replicate within cellular vesicles and transition into non-purulent biofilms. Pneumococci within mature cardiac microlesions exhibited salient biofilm features including intrinsic resistance to antibiotic killing and the presence of an extracellular matrix. Dual RNA-seq and subsequent principal component analyses of heart- and blood-isolated pneumococci confirmed the biofilm phenotype in vivo and revealed stark anatomical site-specific differences in virulence gene expression; the latter having major implications on future vaccine antigen selection. Our RNA-seq approach also identified three genomic islands as exclusively expressed in vivo. Deletion of one such island, Region of Diversity 12, resulted in a biofilm-deficient and highly inflammogenic phenotype within the heart; indicating a possible link between the biofilm phenotype and a dampened host-response. We subsequently determined that biofilm pneumococci released greater amounts of the toxin pneumolysin than did planktonic or RD12 deficient pneumococci. This allowed heart-invaded wildtype pneumococci to kill resident cardiac macrophages and subsequently subvert cytokine/chemokine production and neutrophil infiltration into the myocardium. This is the first report for pneumococcal biofilm formation in an invasive disease setting. We show that biofilm pneumococci actively suppress the host response through pneumolysin-mediated immune cell killing. As such, our findings contradict the emerging notion that biofilm pneumococci are passively immunoquiescent.

Aligner optimization increases accuracy and decreases compute times in multi-species sequence data

As sequencing technologies have evolved, the tools to analyze these sequences have made similar advances. However, for multi-species samples, we observed important and adverse differences in alignment specificity and computation time for bwa- mem (Burrows-Wheeler aligner-maximum exact matches) relative to bwa-aln. Therefore, we sought to optimize bwa-mem for alignment of data from multi-species samples in order to reduce alignment time and increase the specificity of alignments. In the multi-species cases examined, there was one majority member (i.e. Plasmodium falciparum or Brugia malayi) and one minority member (i.e. human or the Wolbachia endosymbiont wBm) of the sequence data. Increasing bwa-mem seed length from the default value reduced the number of read pairs from the majority sequence member that incorrectly aligned to the reference genome of the minority sequence member. Combining both source genomes into a single reference genome increased the specificity of mapping, while also reducing the central processing unit (CPU) time. In Plasmodium, at a seed length of 18 nt, 24.1 % of reads mapped to the human genome using 1.7±0.1 CPU hours, while 83.6 % of reads mapped to the Plasmodium genome using 0.2±0.0 CPU hours (total: 107.7 % reads mapping; in 1.9±0.1 CPU hours). In contrast, 97.1 % of the reads mapped to a combined Plasmodium-human reference in only 0.7±0.0 CPU hours. Overall, the results suggest that combining all references into a single reference database and using a 23 nt seed length reduces the computational time, while maximizing specificity. Similar results were found for simulated sequence reads from a mock metagenomic data set. We found similar improvements to computation time in a publicly available human-only data set.

Effect of a small molecule Lipid II binder on bacterial cell wall stress

We have recently identified small molecule compounds that act as binders of Lipid II, an essential precursor of bacterial cell wall biosynthesis. Lipid II comprised a hydrophilic head group that includes a peptidoglycan subunit composed of N-acetylglucosamine (GlcNAc) and N-acetylmuramic acid (MurNAc) coupled to a short pentapeptide moiety. This headgroup is coupled to a long bactoprenol chain via a pyrophosphate group. Here, we report on the cell wall activity relationship of dimethyl-3-methyl(phenyl)amino-ethenylcyclohexylidene-propenyl-3-ethyl-1,3-benzothiazolium iodide (compound 5107930) obtained by functional and genetic analyses. Our results indicate that compounds bind to Lipid II and cause specific upregulation of the vancomycin-resistance associated gene vraX. vraX is implicated in the cell wall stress stimulon that confers glycopeptide resistance. Our small molecule Lipid II inhibitor retained activity against strains of Staphylococcus aureus mutated in genes encoding the cell wall stress stimulon. This suggests the feasibility of developing this new scaffold as a therapeutic agent in view of increasing glycopeptide resistance.

Genome-wide diversity and gene expression profiling of Babesia microti isolates identify polymorphic genes that mediate host-pathogen interactions

Babesia microti, a tick-transmitted, intraerythrocytic protozoan parasite circulating mainly among small mammals, is the primary cause of human babesiosis. While most cases are transmitted by Ixodes ticks, the disease may also be transmitted through blood transfusion and perinatally. A comprehensive analysis of genome composition, genetic diversity, and gene expression profiling of seven B. microti isolates revealed that genetic variation in isolates from the Northeast United States is almost exclusively associated with genes encoding the surface proteome and secretome of the parasite. Furthermore, we found that polymorphism is restricted to a small number of genes, which are highly expressed during infection. In order to identify pathogen-encoded factors involved in host-parasite interactions, we screened a proteome array comprised of 174 B. microti proteins, including several predicted members of the parasite secretome. Using this immuno-proteomic approach we identified several novel antigens that trigger strong host immune responses during the onset of infection. The genomic and immunological data presented herein provide the first insights into the determinants of B. microti interaction with its mammalian hosts and their relevance for understanding the selective pressures acting on parasite evolution.

Genome-Wide Anaplasma phagocytophilum AnkA-DNA Interactions Are Enriched in Intergenic Regions and Gene Promoters and Correlate with Infection-Induced Differential Gene Expression

Anaplasma phagocytophilum, an obligate intracellular prokaryote, infects neutrophils, and alters cardinal functions via reprogrammed transcription. Large contiguous regions of neutrophil chromosomes are differentially expressed during infection. Secreted A. phagocytophilum effector AnkA transits into the neutrophil or granulocyte nucleus to complex with DNA in heterochromatin across all chromosomes. AnkA binds to gene promoters to dampen cis-transcription and also has features of matrix attachment region (MAR)-binding proteins that regulate three-dimensional chromatin architecture and coordinate transcriptional programs encoded in topologically-associated chromatin domains. We hypothesize that identification of additional AnkA binding sites will better delineate how A. phagocytophilum infection results in reprogramming of the neutrophil genome. Using AnkA-binding ChIP-seq, we showed that AnkA binds broadly throughout all chromosomes in a reproducible pattern, especially at: (i) intergenic regions predicted to be MARs; (ii) within predicted lamina-associated domains; and (iii) at promoters ≤ 3000 bp upstream of transcriptional start sites. These findings provide genome-wide support for AnkA as a regulator of cis-gene transcription. Moreover, the dominant mark of AnkA in distal intergenic regions known to be AT-enriched, coupled with frequent enrichment in the nuclear lamina, provides strong support for its role as a MAR-binding protein and genome “re-organizer.” AnkA must be considered a prime candidate to promote neutrophil reprogramming and subsequent functional changes that belie improved microbial fitness and pathogenicity.

An integrated genomic and transcriptomic survey of mucormycosis-causing fungi

Mucormycosis is a life-threatening infection caused by Mucorales fungi. Here we sequence 30 fungal genomes, and perform transcriptomics with three representative Rhizopus and Mucor strains and with human airway epithelial cells during fungal invasion, to reveal key host and fungal determinants contributing to pathogenesis. Analysis of the host transcriptional response to Mucorales reveals platelet-derived growth factor receptor B (PDGFRB) signaling as part of a core response to divergent pathogenic fungi; inhibition of PDGFRB reduces Mucorales-induced damage to host cells. The unique presence of CotH invasins in all invasive Mucorales, and the correlation between CotH gene copy number and clinical prevalence, are consistent with an important role for these proteins in mucormycosis pathogenesis. Our work provides insight into the evolution of this medically and economically important group of fungi, and identifies several molecular pathways that might be exploited as potential therapeutic targets.

Simultaneous transcriptional profiling of bacteria and their host cells

We developed an RNA-Seq-based method to simultaneously capture prokaryotic and eukaryotic expression profiles of cells infected with intracellular bacteria. As proof of principle, this method was applied to Chlamydia trachomatis-infected epithelial cell monolayers in vitro, successfully obtaining transcriptomes of both C. trachomatis and the host cells at 1 and 24 hours post-infection. Chlamydiae are obligate intracellular bacterial pathogens that cause a range of mammalian diseases. In humans chlamydiae are responsible for the most common sexually transmitted bacterial infections and trachoma (infectious blindness). Disease arises by adverse host inflammatory reactions that induce tissue damage & scarring. However, little is known about the mechanisms underlying these outcomes. Chlamydia are genetically intractable as replication outside of the host cell is not yet possible and there are no practical tools for routine genetic manipulation, making genome-scale approaches critical. The early timeframe of infection is poorly understood and the host transcriptional response to chlamydial infection is not well defined. Our simultaneous RNA-Seq method was applied to a simplified in vitro model of chlamydial infection. We discovered a possible chlamydial strategy for early iron acquisition, putative immune dampening effects of chlamydial infection on the host cell, and present a hypothesis for Chlamydia-induced fibrotic scarring through runaway positive feedback loops. In general, simultaneous RNA-Seq helps to reveal the complex interplay between invading bacterial pathogens and their host mammalian cells and is immediately applicable to any bacteria/host cell interaction.

Disruption of RAB40AL function leads to Martin--Probst syndrome, a rare X-linked multisystem neurodevelopmental human disorder

BACKGROUND AND AIM

Martin--Probst syndrome (MPS) is a rare X-linked disorder characterised by deafness, cognitive impairment, short stature and distinct craniofacial dysmorphisms, among other features. The authors sought to identify the causative mutation for MPS.

METHODS AND RESULTS

Massively parallel sequencing in two affected, related male subjects with MPS identified a RAB40AL (also called RLGP) missense mutation (chrX:102,079,078-102,079,079AC→GA p.D59G; hg18). RAB40AL encodes a small Ras-like GTPase protein with one suppressor of cytokine signalling box. The p.D59G variant is located in a highly conserved region of the GTPase domain between β-2 and β-3 strands. Using RT-PCR, the authors show that RAB40AL is expressed in human fetal and adult brain and kidney, and adult lung, heart, liver and skeletal muscle. RAB40AL appears to be a primate innovation, with no orthologues found in mouse, Xenopus or zebrafish. Western analysis and fluorescence microscopy of GFP-tagged RAB40AL constructs from transiently transfected COS7 cells show that the D59G missense change renders RAB40AL unstable and disrupts its cytoplasmic localisation.

CONCLUSIONS

This is the first study to show that mutation of RAB40AL is associated with a human disorder. Identification of RAB40AL as the gene mutated in MPS allows for further investigations into the molecular mechanism(s) of RAB40AL and its roles in diverse processes such as cognition, hearing and skeletal development.

Identification of rare X-linked neuroligin variants by massively parallel sequencing in males with autism spectrum disorder

BACKGROUND

Autism spectrum disorder (ASD) is highly heritable, but the genetic risk factors for it remain largely unknown. Although structural variants with large effect sizes may explain up to 15% ASD, genome-wide association studies have failed to uncover common single nucleotide variants with large effects on phenotype. The focus within ASD genetics is now shifting to the examination of rare sequence variants of modest effect, which is most often achieved via exome selection and sequencing. This strategy has indeed identified some rare candidate variants; however, the approach does not capture the full spectrum of genetic variation that might contribute to the phenotype.

METHODS

We surveyed two loci with known rare variants that contribute to ASD, the X-linked neuroligin genes by performing massively parallel Illumina sequencing of the coding and noncoding regions from these genes in males from families with multiplex autism. We annotated all variant sites and functionally tested a subset to identify other rare mutations contributing to ASD susceptibility.

RESULTS

We found seven rare variants at evolutionary conserved sites in our study population. Functional analyses of the three 3' UTR variants did not show statistically significant effects on the expression of NLGN3 and NLGN4X. In addition, we identified two NLGN3 intronic variants located within conserved transcription factor binding sites that could potentially affect gene regulation.

CONCLUSIONS

These data demonstrate the power of massively parallel, targeted sequencing studies of affected individuals for identifying rare, potentially disease-contributing variation. However, they also point out the challenges and limitations of current methods of direct functional testing of rare variants and the difficulties of identifying alleles with modest effects.

Microdeletions of 3q29 confer high risk for schizophrenia

Schizophrenia (SZ) is a severe psychiatric illness that affects approximately 1% of the population and has a strong genetic underpinning. Recently, genome-wide analysis of copy-number variation (CNV) has implicated rare and de novo events as important in SZ. Here, we report a genome-wide analysis of 245 SZ cases and 490 controls, all of Ashkenazi Jewish descent. Because many studies have found an excess burden of large, rare deletions in cases, we limited our analysis to deletions over 500 kb in size. We observed seven large, rare deletions in cases, with 57% of these being de novo. We focused on one 836 kb de novo deletion at chromosome 3q29 that falls within a 1.3-1.6 Mb deletion previously identified in children with intellectual disability (ID) and autism, because increasing evidence suggests an overlap of specific rare copy-number variants (CNVs) between autism and SZ. By combining our data with prior CNV studies of SZ and analysis of the data of the Genetic Association Information Network (GAIN), we identified six 3q29 deletions among 7545 schizophrenic subjects and one among 39,748 controls, resulting in a statistically significant association with SZ (p = 0.02) and an odds ratio estimate of 17 (95% confidence interval: 1.36-1198.4). Moreover, this 3q29 deletion region contains two linkage peaks from prior SZ family studies, and the minimal deletion interval implicates 20 annotated genes, including PAK2 and DLG1, both paralogous to X-linked ID genes and now strong candidates for SZ susceptibility.