Faculty Bibliography

Background/Purpose : Although a number of genetic factors, including susceptibility alleles, have been identified in rheumatoid arthritis (RA), the concordance rate in monozygotic (MZ) twins is only 15% (Silman AJ et al, Br J Rheumatol 1993). Given this relatively low rate, environmental factors (including smoking, the microbiome and others) likely play a significant role in disease pathogenesis. Several prior studies have highlighted this relationship. Two examples include oral Porphyromonas gingivalis , which has been implicated in the citrullination of peptides (Hitchon CA et al, J Rheumatol 2010), and intestinal Prevotella copri , which is significantly increased in new-onset RA (Scher JU et al, Elife, 2013). In our study, we sought to further understand this relationship by exploring the microbial, metabolomic and cytokine differences in MZ twins with discordant disease. Methods : Fecal and blood samples were collected from nine pairs of MZ twins where one sibling was diagnosed with RA and the other unaffected. Fecal samples underwent bacterial DNA extraction, amplification and 16S rRNA gene sequencing. Additionally, gas chromatography mass spectrometry (GC-MS) was used to quantify fecal metabolites and multiplex assays were used to quantify fecal and plasma ACPA autoantigens as well as other cytokines (Sokolove J et al, PLoS One, 2012.). Analysis was performed using R, Quantitative Insights into Microbial Ecology (QIIME) and Linear discriminant analysis Effect Size (LEfSe). Results : Microbiome analysis revealed no significant differences in overall bacterial alpha or beta diversity between unaffected and RA twins. On average, RA twins had higher relative abundance of Bacteroidales (RA 56.6% vs Unaffected 47.9%) and lower abundance of Clostridiales (RA 34.5% vs Unaffected 42.8%). LEfSe analysis showed differentially higher abundance of Unclassified Veillonellaceae in unaffected twins (p=0.042). Unaffected twins also demonstrated higher levels of fecal octanoate (p=0.008), a medium-chain fatty acid with beneficial immune effects, as well as significantly increased levels of plasma IL-3 (p=0.039). Conversely, several fecal and plasma citrullinated/ non-citrullinated peptides were significantly higher in RA twins compared to their unaffected siblings (p< 0.05). Conclusion : We characterized differences in the bacterial microbiome, metabolites, cytokines and citrullinated/noncitrullinated peptides in MZ twins discordant for RA. We found changes in the abundance of particular taxa, as well as higher levels of octanoate and IL-3 in unaffected twins, further suggesting the possibility that the initial citrullination and autoimmune events in this disease may occur in the intestinal mucosa. Understanding the immune mechanisms Figure 2 Unaffected twins demonstrate higher levels of octanoate and IL-3, whereas RA twins have significantly higher levels of fecal and plasma citrullinated/non-citrullinated peptides. (A-B) Boxplots of fatty acid metabolites and IL-3 abundance comparing unaffected and RA twins. (C) Tables of plasma and fecal citrullinated/non-citrullinated peptides comparing unaffected and RA twins. Only the top ten peptides by p-value are listed. For all panels, statistical significance calculated using the Wilcoxon signed-rank test. orchestrated by the gut microbiota and their downstream effects may ultimately shed light on the pathogenesis of RA beyond genetic susceptibility and allow us to potentially alter the course of disease development and progression

Lung cancer remains the leading cause of cancer death worldwide1. With new treatment modalities, there has been a shift in focus to how we can predict who may respond to targeted treatments. Current data suggest that the human microbiome can affect lung cancer treatment through its effects on the systemic immune tone. Our group has shown that the lower airway microbiota of lung cancer patients is characterized by enrichment with oral commensals2 which triggers transcriptomic signatures (PI3K, MAPK) previously described in NSCLC 2,3. The impact of local lung dysbiosis on lung cancer progression and survival is unknown. Patients with suspicious nodules on imaging who underwent bronchoscopy were recruited. High-throughput sequencing of bacterial 16S rRNA-encoding gene amplicons was performed. Clustering was based on Dirichlet-Multinomial mixtures (DMM) modeling. RNAseq was performed on bronchial epithelial cells obtained by airway brushing. We focused our analysis on 83 NSCLC samples. Overall alpha-diversity showed that advanced stage (IIIb-VI) lower airway samples were more similar to buccal samples than local stage (I-IIIa), p<0.0001. In addition, worse 6-month and 1-year survival was associated with more similar alpha-diversity between lower airway and buccal samples (Figure 1A-D). Utilizing DMM two clusters were identified, Supraglottic-Predominant-Taxa (SPT) and Background-Predominant-Taxa (BPT). There was a significant increase in percentage of SPT in advance stage compared to local stage (p<0.008) Kaplan-Meir survival analysis shows worse survival in those with NSCLC who were clustered into the SPT group compared to BPT (p=0.0003, Figure 1E). With RNAseq, differentially expressed genes between advanced stage vs. local stage and 6-month vs. 1-year survival were not as pronounced as SPT vs. BPT (Figure 1F) suggesting that globally, transcriptomic changes between different stage and NSCLC survival were difficult to detect as compared to when airway microbiome were differentiated. In lung cancer, dysbiosis within the lower airway microenvironment, possibly by micro-aspiration, is associated with a worse 6-month and 1-year survival. This change is also associated with transcriptome changes in the local environment

BACKGROUND: While most hemangioblastomas (~70%) are sporadic and occur predominantly in the cerebellum, they may present as well as familial form associated with von Hippel-Lindau (VHL) syndrome, an autosomal dominant disorder caused by germline mutations of the VHL gene that trigger nuclear translocation of hypoxia-inducible factor (HIF)- 1alpha and angiogenesis. Although inactivation of VHL, a tumor suppressor gene, has been observed in hemangioblastomas, the underlying pathogenic mechanisms responsible for familial and sporadic hemangioblastomas remain incompletely understood.
METHOD(S): Whole exome sequencing of cerebellar hemangioblastoma tumors and matched blood leukocytes from 24 patients, age 24-63, was performed. After preparation and amplification of barcoded libraries, exomes were captured using Kapa Biosystems methodology and paired-end sequenced on Illumina HiSeq 2500 to an average 100-fold coverage. Following read alignment to hg19 genome, tumor and germline (leukocyte) sequences were compared, and pathogenic single nucleotide variants (SNVs) identified and validated by re-sequencing followed by pathway analysis. Additionally, tumor RNA isolated using Maxwell Promega was sequenced on Illumina instrument and the expression counts determined and normalized.
RESULT(S): We found 314 pathogenic and/or highly deleterious mutations (both germline and somatic) with a median of 13 mutations per patient. Five patients had VHL syndrome (germline VHL mutation) and 4 carried somatic VHL mutations. Among the VHL tumors, 82 mutations were identified, including HNF1B, NOTCH1 and TCF7L1, suggesting a potential contribution of altered RNA metabolism based upon pathway analysis. Among all hemangioblastomas, germline growth factor receptor variants (FGFR4 p.G388R (14/23 (61%) patients), IGF1R, PDGFRA and TYK2) known to activate STAT3 signaling and induce HIF-1alpha and angiogenesis, were identified. Non-hierarchical clustering of RNA sequencing data revealed two transcriptionally-distinct subtypes of hemangioblastomas.
CONCLUSION(S): Our findings indicate that hemangioblastomas can also occur by germline mutations known to activate STAT3 signaling, which may have significant implication in genetic testing and counseling of patients with hemangioblastomas

Background: Plakophilin-2 (PKP2) is classically defined as a protein of the desmosome, an intercellular adhesion structure that also acts as a signaling hub to maintain structural and electrical homeostasis. Mutations in PKP2 associate with most cases of gene-positive arrhythmogenic right ventricular cardiomyopathy (ARVC). A better understanding of PKP2 cardiac biology can help elucidate the mechanisms underlying arrhythmic and cardiomyopathic events that occur consequent to its mutation. Here we sought to captureearly molecular/cellular events that can act as nascent substrates for subsequent arrhythmic/cardiomyopathic phenotypes.
Method(s): We used multiple quantitative imaging modalities, as well as biochemical and high-resolution mass spectrometry methods to study the functional/structural properties of cells/tissues derived from cardiomyocytespecific, tamoxifen-activated, PKP2 knockout mice ("PKP2cKO"). Studies were carried out 14 days post-tamoxifen injection, a time point preceding an overt electrical or structural phenotype.Myocytes from right or left ventricular free wall were studied separately, to detect functional/structural asymmetries.
Result(s): Most properties of PKP2cKO left ventricular (LV) myocytes were not different from control; in contrast, PKP2cKO right ventricular (RV) myocytes showed increased amplitude and duration of Ca2+transients, increased frequency of spontaneous Ca2+release events, increased [Ca2+] in the cytoplasm and sarcoplasmic reticulum compartments, and dynamic Ca2+accumulation in mitochondria. In addition, RyR2 in RV presented enhanced sensitivity to Ca2+and preferential phosphorylation in a domain known to modulate Ca2+gating. RNAseq at 14 days post-TAM showed no relevant difference in transcript abundance between RV and LV, neither in control nor in PKP2cKO cells, suggesting that in the earliest stage, [Ca2+]i dysfunction is not transcriptional. Rather, we found an RV-predominant increase in membrane permeability that can permit Ca2+entry into the cell. Cx43 ablation mitigated the increase in membrane permeability, the accumulation of cytoplasmic Ca2+and the early stages of RV dysfunction.
Conclusion(s): Loss of PKP2 creates an RV-predominant arrhythmogenic substrate (Ca2+ dysregulation) that precedes the cardiomyopathy and that is, at least in part, mediated by a Cx43-dependent membrane conduit. Given that asymmetric Ca2+ dysregulation precedes the cardiomyopathic stage, we speculate that abnormal Ca2+ handling in RV myocytes can be a trigger for gross structural changes observed at a later stage

Local translation can support memory consolidation by supplying new proteins to synapses undergoing plasticity. Translation in adult forebrain dendrites is an established mechanism of synaptic plasticity and is regulated by learning, yet there is no evidence for learning-regulated protein synthesis in adult forebrain axons, which have traditionally been believed to be incapable of translation. Here we show that axons in the adult rat amygdala contain translation machinery, and use translating ribosome affinity purification (TRAP) with RNASeq to identify mRNAs in cortical axons projecting to the amygdala, over 1200 of which were regulated during consolidation of associative memory. Mitochondrial and translation-related genes were upregulated, whereas synaptic, cytoskeletal, and myelin-related genes were downregulated; the opposite effects were observed in the cortex. Our results demonstrate that axonal translation occurs in the adult forebrain and is altered after learning, supporting the likelihood that local translation is more a rule than an exception in neuronal processes.

Multifocal breast cancer (MFBC), ductal type, has been hypothesized to arise by one of two mechanisms: either through intramammary/intralymphatic spread from a single index tumor (MBC-1), or as multiple independent tumors with each focus carrying its corresponding ductal carcinoma in-situ (MBC-2). In order to improve our understanding of MFBC pathogenesis, we employed laser capture microdissection coupled with whole-exome sequencing to study clonal origin in MFBC. We selected three cases of MBC-1 (C1 to C3) and MBC-2 (C4 to C6) and analyzed three foci from each case. MBC-1 cases were histologically similar and showed a strong predilection for satellite foci, vascular invasion and nodal metastasis when compared to MBC-2. Our bioinformatics approach provided strong evidence for clonal relationships in MBC-1, as demonstrated by distinct clusters of genes conserved across all tumor foci. Conversely, no gene clusters were shared across all the foci in MBC-2, suggesting multiple independent tumors. These findings provide further support for the two distinct pathogenetic mechanisms in MFBC.

Recurrent somatic missense mutations in histone H3 genes have been identified in subsets of pediatric cancers. H3K36 histone mutations have recently been recognized as oncogenic drivers in rare subsets of malignant soft tissue sarcomas but have not been reported in histiocytic neoplasms. Currently, the histological and molecular spectrum, as well as the clinical behavior of H3K36-mutant soft tissue malignancies, is largely unknown. We describe a pediatric patient with a HIST1H3B K36I-mutant histiocytic tumor arising in the skull. After the failure of upfront therapy for histiocytosis and development of widely disseminated metastatic disease, the patient had an exceptional response to empiric chemotherapy and remains in complete disease remission for more than 5 years. Our report expands the histological spectrum of H3K36M/I-mutant soft tissue malignancies to histiocytic neoplasms and indicates that multiagent sarcoma-like chemotherapy can be highly effective even in the setting of widely disseminated metastatic disease.

[This corrects the article DOI: 10.1371/journal.pone.0212593.].