Loss of developmentally derived Irf8+ macrophages promotes hyperinnervation and arrhythmia in the adult zebrafish heart

Recent developments in cardiac macrophage biology have broadened our understanding of the critical functions of macrophages in the heart. As a result, there is further interest in understanding the independent contributions of distinct subsets of macrophage to cardiac development and function. Here, we demonstrate that genetic loss of interferon regulatory factor 8 (Irf8)-positive embryonic-derived macrophages significantly disrupts cardiac conduction, chamber function, and innervation in adult zebrafish. At 4 months post-fertilization (mpf), homozygous irf8st96/st96 mutants have significantly shortened atrial action potential duration and significant differential expression of genes involved in cardiac contraction. Functional in vivo assessments via electro- and echocardiograms at 12 mpf reveal that irf8 mutants are arrhythmogenic and exhibit diastolic dysfunction and ventricular stiffening. To identify the molecular drivers of the functional disturbances in irf8 null zebrafish, we perform single cell RNA sequencing and immunohistochemistry, which reveal increased leukocyte infiltration, epicardial activation, mesenchymal gene expression, and fibrosis. Irf8 null hearts are also hyperinnervated and have aberrant axonal patterning, a phenotype not previously assessed in the context of cardiac macrophage loss. Gene ontology analysis supports a novel role for activated epicardial-derived cells (EPDCs) in promoting neurogenesis and neuronal remodeling in vivo. Together, these data uncover significant cardiac abnormalities following embryonic macrophage loss and expand our knowledge of critical macrophage functions in heart physiology and governing homeostatic heart health.

Single cell transcriptional changes across the blood stages of artemisinin resistant K13C580Y mutant Plasmodium falciparum upon dihydroartemisinin exposure

Artemisinins have been a cornerstone of malaria control, but resistance in Plasmodium falciparum, due to mutations in the Kelch 13 gene, threaten these advances. Artemisinin exposure results in a dynamic transcriptional response across multiple pathways, but most work has focused on ring stages and ex vivo transcriptional analysis, limiting evaluation of all life cycle stages. We applied single cell RNAseq to two unsynchronized isogenic parasite lines (K13C580 and K13580Y) over 6 hrs after a pulse exposure to dihydroartemisinin (DHA). Transcription was altered across all stages, with the greatest occurring at the early trophozoite and mid ring stage in both lines. This response involved the arrest of metabolic processes and the enhancement of protein trafficking and the unfolded protein response. While similar, the response was enhanced in the K13580Y mutant, which may lead to the dormancy phenomenon upon treatment. Increased surface protein expression was seen in mutant parasites at baseline and upon drug exposure, highlighted by the increased expression of PfEMP1 and GARP, a potential therapeutic target. Antibody targeting GARP maintained anti-parasitic efficacy in mutant parasites. This work provides single cell insight of gene transcription across all life cycle stages revealing transcriptional changes that could initiate dormancy state and mediate survival.

Viral and host factors drive a type 1 Epstein-Barr virus spontaneous lytic phenotype

IMPORTANCE

Epstein-Barr virus (EBV) infects over 95% of adults worldwide. Given its connection to various cancers and autoimmune disorders, it is important to understand the mechanisms by which infection with EBV can lead to these diseases. In this study, we describe an unusual spontaneous lytic phenotype in EBV strains isolated from Kenyan endemic Burkitt lymphoma patients. Because lytic replication of EBV has been linked to the pathogenesis of various diseases, these data could illuminate viral and host factors involved in this process.

Expression Microdissection for the Analysis of miRNA in a Single-Cell Type

Cell-specific microRNA (miRNA) expression estimates are important in characterizing the localization of miRNA signaling within tissues. Much of these data are obtained from cultured cells, a process known to significantly alter miRNA expression levels. Thus, our knowledge of in vivo cell miRNA expression estimates is poor. We previously demonstrated expression microdissection-miRNA-sequencing (xMD-miRNA-seq) to acquire in vivo estimates, directly from formalin-fixed tissues, albeit with a limited yield. In this study, we optimized each step of the xMD process, including tissue retrieval, tissue transfer, film preparation, and RNA isolation, to increase RNA yields and ultimately show strong enrichment for in vivo miRNA expression by qPCR array. These method improvements, such as the development of a noncrosslinked ethylene vinyl acetate membrane, resulted in a 23- to 45-fold increase in miRNA yield, depending on the cell type. By qPCR, miR-200a increased by 14-fold in xMD-derived small intestine epithelial cells, with a concurrent 336-fold reduction in miR-143 relative to the matched nondissected duodenal tissue. xMD is now an optimized method to obtain robust in vivo miRNA expression estimates from cells. xMD will allow formalin-fixed tissues from surgical pathology archives to make theragnostic biomarker discoveries.

Endemic Burkitt lymphoma avatar mouse models for exploring inter-patient tumor variation and testing targeted therapies

Endemic Burkitt lymphoma (BL) is a childhood cancer in sub-Saharan Africa characterized by Epstein-Barr virus and malaria-associated aberrant B-cell activation and MYC chromosomal translocation. Survival rates hover at 50% after conventional chemotherapies; therefore, clinically relevant models are necessary to test additional therapies. Hence, we established five patient-derived BL tumor cell lines and corresponding NSG-BL avatar mouse models. Transcriptomics confirmed that our BL lines maintained fidelity from patient tumors to NSG-BL tumors. However, we found significant variation in tumor growth and survival among NSG-BL avatars and in Epstein-Barr virus protein expression patterns. We tested rituximab responsiveness and found one NSG-BL model exhibiting direct sensitivity, characterized by apoptotic gene expression counterbalanced by unfolded protein response and mTOR pro-survival pathways. In rituximab-unresponsive tumors, we observed an IFN-α signature confirmed by the expression of IRF7 and ISG15. Our results demonstrate significant inter-patient tumor variation and heterogeneity, and that contemporary patient-derived BL cell lines and NSG-BL avatars are feasible tools to guide new therapeutic strategies and improve outcomes for these children.

Distinctive Kaposi Sarcoma-Associated Herpesvirus Serological Profile during Acute Plasmodium falciparum Malaria Episodes

The seroprevalence of Kaposi sarcoma-associated herpesvirus (KSHV) and the incidence of endemic Kaposi sarcoma (KS) overlap with regions of malaria endemicity in sub-Saharan Africa. Multiple studies have shown an increased risk of KSHV seroconversion in children from high malaria compared to low malaria regions; however, the impact of acute episodes of Plasmodium falciparum (P. falciparum) malaria on KSHV's biphasic life cycle and lytic reactivation has not been determined. Here, we examined KSHV serological profiles and viral loads in 134 children with acute malaria and 221 healthy children from high malaria regions in Kisumu, as well as 77 healthy children from low malaria regions in Nandi. We assayed KSHV, Epstein-Barr virus (EBV), and P. falciparum malaria antibody responses in these three by multiplexed Luminex assay. We confirmed that KSHV seroprevalence was significantly associated with malaria endemicity (OR = 1.95, 1.18-3.24 95% CI, p = 0.01) with 71-77% seropositivity in high-malaria (Kisumu) compared to 28% in low-malaria (Nandi) regions. Furthermore, KSHV serological profiles during acute malaria episodes were distinct from age-matched non-malaria-infected children from the same region. Paired IgG levels also varied after malaria treatment, with significantly higher anti-ORF59 at day 0 but elevated ORF38, ORF73, and K8.1 at day 3. Acute malaria episodes is characterized by perturbation of KSHV latency in seropositive children, providing further evidence that malaria endemicity contributes to the observed increase in endemic KS incidence in sub-Saharan Africa.

Profiling genome-wide recombination in Epstein Barr virus reveals type-specific patterns and associations with endemic-Burkitt lymphoma

BACKGROUND

Endemic Burkitt lymphoma (eBL) is potentiated through the interplay of Epstein Barr virus (EBV) and holoendemic Plasmodium falciparum malaria. To better understand EBV's biology and role in eBL, we characterized genome-wide recombination sites and patterns as a source of genetic diversity in EBV genomes in our well-defined population of eBL cases and controls from Western Kenya.

METHODS

EBV genomes representing 54 eBL cases and 32 healthy children from the same geographic region in Western Kenya that we previously sequenced were analyzed. Whole-genome multiple sequence alignment, recombination analyses, and phylogenetic inference were made using multiple alignment with fast Fourier transform, recombination detection program 4, and molecular evolutionary genetics analysis.

RESULTS

We identified 28 different recombination events and 71 (82.6%) of the 86 EBV genomes analyzed contained evidence of one or more recombinant segments. Associated recombination breakpoints were found to occur in a total of 42 different genes, with only 7 (16.67%) being latent genes. Recombination events were major drivers of clustering within genome-wide phylogenetic trees. The occurrence of recombination segments was comparable between genomes from male and female participants and across age groups. More recombinant segments were found in EBV type 1 genomes (p = 6.4e - 06) and the genomes from the eBLs (p = 0.037). Two recombination events were enriched in the eBLs; event 47 (OR = 4.07, p = 0.038) and event 50 (OR = 14.24, p = 0.012).

CONCLUSIONS

EBV genomes have extensive evidence of recombination likely acquired progressively and cumulatively over time. Recombination patterns display a heterogeneous occurrence rate across the genome with enrichment in lytic genes. Overall, recombination appears to be a major evolutionary force impacting EBV diversity and genome structure with evidence of the association of specific recombinants with eBL.

Kaposi Sarcoma-Associated Herpesvirus Infection and Endemic Burkitt Lymphoma

BACKGROUND

Endemic Burkitt lymphoma (eBL) is associated with Epstein-Barr virus (EBV) and Plasmodium falciparum malaria coinfections. However, the role of Kaposi sarcoma-associated herpesvirus (KSHV), also endemic in Africa, has not been evaluated as a cofactor in eBL pathogenesis.

METHODS

Multiplexed seroprofiles for EBV, malaria, and KSHV were generated for 266 eBL patients, 78 non-eBL cancers, and 202 healthy children. KSHV and EBV loads were quantified by PCR.

RESULTS

KSHV seroprevalence did not differ by study group but was associated with age. Seropositivity, defined by K8.1/LANA or in combination with 5 other KSHV antigens (ORF59, ORF65, ORF61, ORF38, and K5) was associated with antimalarial antibody levels to AMA1 (odds ratio [OR], 2.41, P < .001; OR, 2.07, P < .001) and MSP1 (OR, 2.41, P = .0006; OR, 5.78, P < .001), respectively. KSHV loads did not correlate with antibody levels nor differ across groups but were significantly lower in children with detectable EBV viremia (P = .014).

CONCLUSIONS

Although KSHV-EBV dual infection does not increase eBL risk, EBV appears to suppress reactivation of KSHV while malaria exposure is associated with KSHV infection and/or reactivation. Both EBV and malaria should, therefore, be considered as potential effect modifiers for KSHV-associated cancers in sub-Saharan Africa.

A New Hope for CD56negCD16pos NK Cells as Unconventional Cytotoxic Mediators: An Adaptation to Chronic Diseases

Natural Killer (NK) cells play an essential role in antiviral and anti-tumoral immune responses. In peripheral blood, NK cells are commonly classified into two major subsets: CD56^brightCD16^neg and CD56^dimCD16^pos despite the characterization of a CD56^negCD16^pos subset 25 years ago. Since then, several studies have described the prevalence of an CD56^negCD16^pos NK cell subset in viral non-controllers as the basis for their NK cell dysfunction. However, the mechanistic basis for their cytotoxic impairment is unclear. Recently, using a strict flow cytometry gating strategy to exclude monocytes, we reported an accumulation of CD56^negCD16^pos NK cells in Plasmodium falciparum malaria-exposed children and pediatric cancer patients diagnosed with endemic Burkitt lymphoma (eBL). Here, we use live-sorted cells, histological staining, bulk RNA-sequencing and flow cytometry to confirm that this CD56^negCD16^pos NK cell subset has the same morphological features as the other NK cell subsets and a similar transcriptional profile compared to CD56^dimCD16^pos NK cells with only 120 genes differentially expressed (fold change of 1.5, p < 0.01 and FDR<0.05) out of 9235 transcripts. CD56^negCD16^pos NK cells have a distinct profile with significantly higher expression of MPEG1 (perforin 2), FCGR3B (CD16b), FCGR2A, and FCGR2B (CD32A and B) as well as CD6, CD84, HLA-DR, LILRB1/2, and PDCD1 (PD-1), whereas Interleukin 18 (IL18) receptor genes (IL18RAP and IL18R1), cytotoxic genes such as KLRF1 (NKp80) and NCR1 (NKp46), and inhibitory HAVCR2 (TIM-3) are significantly down-regulated compared to CD56^dimCD16^pos NK cells. Together, these data confirm that CD56^negCD16^pos cells are legitimate NK cells, yet their transcriptional and protein expression profiles suggest their cytotoxic potential is mediated by pathways reliant on antibodies such as antibody-dependent cell cytotoxicity (ADCC), antibody-dependent respiratory burst (ADRB), and enhanced by complement receptor 3 (CR3) and FAS/FASL interaction. Our findings support the premise that chronic diseases induce NK cell modifications that circumvent proinflammatory mediators involved in direct cytotoxicity. Therefore, individuals with such altered NK cell profiles may respond differently to NK-mediated immunotherapies, infections or vaccines depending on which cytotoxic mechanisms are being engaged.

The whole-genome landscape of Burkitt lymphoma subtypes

Burkitt lymphoma (BL) is an aggressive, MYC-driven lymphoma comprising 3 distinct clinical subtypes: sporadic BLs that occur worldwide, endemic BLs that occur predominantly in sub-Saharan Africa, and immunodeficiency-associated BLs that occur primarily in the setting of HIV. In this study, we comprehensively delineated the genomic basis of BL through whole-genome sequencing (WGS) of 101 tumors representing all 3 subtypes of BL to identify 72 driver genes. These data were additionally informed by CRISPR screens in BL cell lines to functionally annotate the role of oncogenic drivers. Nearly every driver gene was found to have both coding and non-coding mutations, highlighting the importance of WGS for identifying driver events. Our data implicate coding and non-coding mutations in IGLL5, BACH2, SIN3A, and DNMT1. Epstein-Barr virus (EBV) infection was associated with higher mutation load, with type 1 EBV showing a higher mutational burden than type 2 EBV. Although sporadic and immunodeficiency-associated BLs had similar genetic profiles, endemic BLs manifested more frequent mutations in BCL7A and BCL6 and fewer genetic alterations in DNMT1, SNTB2, and CTCF. Silencing mutations in ID3 were a common feature of all 3 subtypes of BL. In vitro, mass spectrometry-based proteomics demonstrated that the ID3 protein binds primarily to TCF3 and TCF4. In vivo knockout of ID3 potentiated the effects of MYC, leading to rapid tumorigenesis and tumor phenotypes consistent with those observed in the human disease.

Comprehensive Transcriptome and Mutational Profiling of Endemic Burkitt Lymphoma Reveals EBV Type-Specific Differences

Endemic Burkitt lymphoma (eBL) is the most common pediatric cancer in malaria-endemic equatorial Africa and nearly always contains Epstein-Barr virus (EBV), unlike sporadic Burkitt lymphoma (sBL) that occurs with a lower incidence in developed countries. Given these differences and the variable clinical presentation and outcomes, we sought to further understand pathogenesis by investigating transcriptomes using RNA sequencing (RNAseq) from multiple primary eBL tumors compared with sBL tumors. Within eBL tumors, minimal expression differences were found based on: anatomical presentation site, in-hospital survival rates, and EBV genome type, suggesting that eBL tumors are homogeneous without marked subtypes. The outstanding difference detected using surrogate variable analysis was the significantly decreased expression of key genes in the immunoproteasome complex (PSMB9/β1i, PSMB10/β2i, PSMB8/β5i, and PSME2/PA28β) in eBL tumors carrying type 2 EBV compared with type 1 EBV. Second, in comparison with previously published pediatric sBL specimens, the majority of the expression and pathway differences was related to the PTEN/PI3K/mTOR signaling pathway and was correlated most strongly with EBV status rather than geographic designation. Third, common mutations were observed significantly less frequently in eBL tumors harboring EBV type 1, with mutation frequencies similar between tumors with EBV type 2 and without EBV. In addition to the previously reported genes, a set of new genes mutated in BL, including TFAP4, MSH6, PRRC2C, BCL7A, FOXO1, PLCG2, PRKDC, RAD50, and RPRD2, were identified. Overall, these data establish that EBV, particularly EBV type 1, supports BL oncogenesis, alleviating the need for certain driver mutations in the human genome.

IMPLICATIONS

Genomic and mutational analyses of Burkitt lymphoma tumors identify key differences based on viral content and clinical outcomes suggesting new avenues for the development of prognostic molecular biomarkers and therapeutic interventions.

Integrative microRNA and mRNA deep-sequencing expression profiling in endemic Burkitt lymphoma

BACKGROUND

Burkitt lymphoma (BL) is characterized by overexpression of the c-myc oncogene, which in the vast majority of cases is a consequence of an IGH/MYC translocation. While myc is the seminal event, BL is a complex amalgam of genetic and epigenetic changes causing dysregulation of both coding and non-coding transcripts. Emerging evidence suggest that abnormal modulation of mRNA transcription via miRNAs might be a significant factor in lymphomagenesis. However, the alterations in these miRNAs and their correlations to their putative mRNA targets have not been extensively studied relative to normal germinal center (GC) B cells.

METHODS

Using more sensitive and specific transcriptome deep sequencing, we compared previously published small miRNA and long mRNA of a set of GC B cells and eBL tumors. MiRWalk2.0 was used to identify the validated target genes for the deregulated miRNAs, which would be important for understanding the regulatory networks associated with eBL development.

RESULTS

We found 211 differentially expressed (DE) genes (79 upregulated and 132 downregulated) and 49 DE miRNAs (22 up-regulated and 27 down-regulated). Gene Set enrichment analysis identified the enrichment of a set of MYC regulated genes. Network propagation-based method and correlated miRNA-mRNA expression analysis identified dysregulated miRNAs, including miR-17~95 cluster members and their target genes, which have diverse oncogenic properties to be critical to eBL lymphomagenesis. Central to all these findings, we observed the downregulation of ATM and NLK genes, which represent important regulators in response to DNA damage in eBL tumor cells. These tumor suppressors were targeted by multiple upregulated miRNAs (miR-19b-3p, miR-26a-5p, miR-30b-5p, miR-92a-5p and miR-27b-3p) which could account for their aberrant expression in eBL.

CONCLUSION

Combined loss of p53 induction and function due to miRNA-mediated regulation of ATM and NLK, together with the upregulation of TFAP4, may be a central role for human miRNAs in eBL oncogenesis. This facilitates survival of eBL tumor cells with the IGH/MYC chromosomal translocation and promotes MYC-induced cell cycle progression, initiating eBL lymphomagenesis. This characterization of miRNA-mRNA interactions in eBL relative to GC B cells provides new insights on miRNA-mediated transcript regulation in eBL, which are potentially useful for new improved therapeutic strategies.

Comprehensive Transcriptome and Mutational Profiling of Endemic Burkitt Lymphoma Reveals EBV Type-Specific Differences

Endemic Burkitt lymphoma (eBL) is the most common pediatric cancer in malaria-endemic equatorial Africa and nearly always contains Epstein-Barr virus (EBV), unlike sporadic Burkitt lymphoma (sBL) that occurs with a lower incidence in developed countries. Given these differences and the variable clinical presentation and outcomes, we sought to further understand pathogenesis by investigating transcriptomes using RNA sequencing (RNAseq) from multiple primary eBL tumors compared with sBL tumors. Within eBL tumors, minimal expression differences were found based on: anatomical presentation site, in-hospital survival rates, and EBV genome type, suggesting that eBL tumors are homogeneous without marked subtypes. The outstanding difference detected using surrogate variable analysis was the significantly decreased expression of key genes in the immunoproteasome complex (PSMB9/β1i, PSMB10/β2i, PSMB8/β5i, and PSME2/PA28β) in eBL tumors carrying type 2 EBV compared with type 1 EBV. Second, in comparison with previously published pediatric sBL specimens, the majority of the expression and pathway differences was related to the PTEN/PI3K/mTOR signaling pathway and was correlated most strongly with EBV status rather than geographic designation. Third, common mutations were observed significantly less frequently in eBL tumors harboring EBV type 1, with mutation frequencies similar between tumors with EBV type 2 and without EBV. In addition to the previously reported genes, a set of new genes mutated in BL, including TFAP4, MSH6, PRRC2C, BCL7A, FOXO1, PLCG2, PRKDC, RAD50, and RPRD2, were identified. Overall, these data establish that EBV, particularly EBV type 1, supports BL oncogenesis, alleviating the need for certain driver mutations in the human genome.

IMPLICATIONS

Genomic and mutational analyses of Burkitt lymphoma tumors identify key differences based on viral content and clinical outcomes suggesting new avenues for the development of prognostic molecular biomarkers and therapeutic interventions.

Human and Epstein-Barr Virus miRNA Profiling as Predictive Biomarkers for Endemic Burkitt Lymphoma

Endemic Burkitt lymphoma (eBL) is an aggressive B cell lymphoma and is associated with Epstein-Barr virus (EBV) and Plasmodium falciparum malaria co-infections. Central to BL oncogenesis is the over-expression of the MYC proto-oncogene which is caused by a translocation of an Ig enhancer in approximation to the myc gene. While whole genome/transcriptome sequencing methods have been used to define driver mutations and transcriptional dysregulation, microRNA (miRNA) dysregulation and differential expression has yet to be fully characterized. We hypothesized that both human and EBV miRNAs contribute to eBL clinical presentation, disease progression, and poor outcomes. Using sensitive and precise deep sequencing, we identified miRNAs from 17 Kenyan eBL patient tumor samples and delineated the complement of both host and EBV miRNAs. One human miRNA, hsa-miR-10a-5p was found to be differentially expressed (DE), being down-regulated in jaw tumors relative to abdominal and in non-survivors compared to survivors. We also examined EBV miRNAs, which made up 2.7% of the miRNA composition in the eBL samples. However, we did not find any significant associations regarding initial patient outcome or anatomical presentation. Gene ontology analysis and pathway enrichment of previously validated targets of miR-10a-5p suggest that it can promote tumor cell survival as well as aid in evasion of apoptosis. To examine miR-10a-5p regulatory effect on gene expression in eBL, we performed a pairwise correlation coefficient analysis on the expression levels of all its validated targets. We found a significant enrichment of correlated target genes consistent with miR-10a-5p impacting expression. The functions of genes and their correlation fit with multiple target genes impacting tumor resilience. The observed downregulation of miR-10a and associated genes suggests a role for miRNA in eBL patient outcomes and has potential as a predictive biomarker that warrants further investigation.

Interleukin-6 and interleukin-10 gene promoter polymorphisms and risk of endemic Burkitt lymphoma

Overexpression of interleukin-6 (IL-6) and IL-10 in endemic Burkitt lymphoma (eBL) may facilitate tumorigenesis by providing a permissive cytokine milieu. Promoter polymorphisms influence interindividual differences in cytokine production. We hypothesized that children genetically predisposed for elevated cytokine levels may be more susceptible to eBL. Using case-control samples from western Kenya consisting of 117 eBL cases and 88 ethnically matched healthy controls, we tested for the association between eBL risk and IL-10 (rs1800896, rs1800871, and rs1800872) and IL-6 (rs1800795) promoter single nucleotide polymorphisms (SNPs) as well as IL-10 promoter haplotypes. In addition, the association between these variants and Epstein Barr Virus (EBV) load was examined. Results showed that selected IL-10 and IL-6 promoter SNPs and IL-10 promoter haplotypes were not associated with risk eBL or EBV levels in EBV-seropositive children. Findings from this study reveal that common variants within the IL-10 and IL-6 promoters do not independently increase eBL risk in this vulnerable population.